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Abstract 

We present a detailed study, in the semi-classical regime /i ^ 0, of microlocal properties 
of systems of two commuting /i-pseudo-differential operators Pi(h), P2{h) such that the joint 
principal symbol p = (pi,P2) has a special kind of singularity called a, focus-focus singularity. 
Typical examples include the quantum spherical pendulum or the quantum Champagne 
bottle. 

In the spirit of Colin de Verdiere and Parisse [0 [l2[ we show that such systems 
have a universal behavior described by singular quantization conditions of Bohr-Sommerfeld 
type. 

These conditions are used to give a precise description of the joint spectrum of such 
systems, including the phenomenon of quantum monodromy and different formulations of 
the counting function for the joint eigenvalues close to the singularity, in which a logarithm of 
the semi-classical constant h appears. Thanks to numerical computations done by M.S. Child 
for the case of the Champagne bottle, we are able to accurately illustrate our statements. 
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1 Introduction 



In the long history of completely integrable systems, an important object was discovered quite 
recently (Duistermaat [^]) : the monodromy of the system, whose non-triviality prevents the 
construction of global action variables. The question about the impact of this invariant on the 



spectrum of quantum integrable system was raised by Cushman and Duistermaat in |15]; an 
answer is proposed in |^^. The issue is to describe the joint spectrum of two commuting h- 
pseudo-differential operators Pi, P2 in a region close to a critical point of the C°° joint principal 
symbol 

when the underlying Liouville integrable system pi , p2 has non trivial monodromy. 

However, because of a well-known drawback of both the usual WKB construction and the 
standard Bohr-Sommerfeld quantization conditions, the descriptions had to keep a reasonable 
distance away from the singular value of p. Here reasonable may be small but means fixed, as 
h tends to zero. In this way, an increasing number of eigenvalues (as /i — > 0) remained out of 
control. 

On the other hand, recent achievements in semi-classical analysis of Schrodinger operators 
near a critical point of the potential, often via the use of microlocal normal forms, suggested that 
this problem should be solvable. I am referring here for instance to the work of Helffer-Robert 



pq| , Helffer-Sjostrand [^], Sjostrand [38|, Marz |33], Brummelhuis-Paul-Uribe y], ... and in 
particular to the articles by Colin de Verdiere and Parisse |11, 13|, in which the case of a 
local maximum of a smooth (C°°) potential for a 1-dimensional Schrodinger operator is treated. 
Their method rested upon a smooth normal form theorem and on the study of the hyperbolicity 
of the classical Hamiltonian flow. 

It turns out that any two degree of freedom quantum integrable system with a focus-focus 
singularity exhibits at the same time a non-trivial monodromy and a hyperbolic behavior of the 
Hamiltonian flow. Moreover, such a singularity admits a smooth normal form, due to Eliasson 



22], that has a semi-classical analogue [^^. This allows us to settle an analysis combining a 
geometrical description of the underlying classical completely integrable system with a microlocal 
analysis near the singularity. 

The main result of this paper is the statement of the singular Bohr-Sommerfeld quantization 
conditions (theorem |6l9| ), which are uniform in a neighborhood of the critical value. In other 
words, these conditions are able to uniformly describe an increasing number of joint eigenvalues, 
and as such they contain the description of the quantum monodromy far from the singularity 
as well as the asymptotic distribution of the joint eigenvalues near the singularity. 

The organization of this paper is as follows. In the first parts, the aim is to settle the 
ingredients needed to have a global picture of the problem. The necessary background concerning 
Liouville integrable systems is recalled (section |2|), and the general notion of semi- classical 
integrable systems is explained in section ^, where we point out the role of both the principal 
and sub-principal symbols. The version of the microlocal tools that we will be using is given 
in section ^; its application to completely integrable systems is developed in the next section 
(H). Besides known results concerning the WKB construction, the most fundamental results 



explaining our methods are probably propositions 5.1 and 5.5. 
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The rest of the paper is devoted to the case of a focus-focus singularity in a 4-dimensional 
cotangent bundle. Section ^ contains the microlocal analysis and the geometrical description 
of the monodromy that lead to our main result (theorem |6.9| ). Finally, we derive in section ^ 
the structure of the joint spectrum. Theorem [7.3| shows that one can count the number of joint 
eigenvalues without knowing their precise distribution near the critical value, and theorem |7.6| , 
on the contrary, studies this distribution, by estimating the gaps between the joint eigenvalues. 
These statements are illustrated by the example of the Champagne bottle. 
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2 Classical completely integrable systems 



2.1 Definition 



Let (M, Lu) be a symplectic manifold of dimension 2n. A completely (or "Liouville" ) integrable 
system on M is the data of n functions /i , . . . , /„ in involution with respect to the symplectic 
Poisson bracket, with the requirement that their differentials dfi are almost everywhere inde- 
pendent. The function 

F :M9m^(/i(m),...,/„(m))GM" 

is called the momentum map. It is indeed a momentum map for the Abelian infinitesimal action 
of R" into the Lie algebra X{M) of Hamiltonian vector fields on M, given by the generators 
{Xi, . . . ,Xn)- Here we have denoted by Xi the Hamiltonian vector field Xf. associated to /j. 
The flows of these vector fields yield a local Abelian action of on M, simply referred to in 
the sequel by "the action" or "the flow" of the system. 

Note that each /j, and hence every function of the form g{fi, ■ ■ ■ , fn) is locally constant 
under this action. The momentum map F thus defines a singular fibration by invariant leaves. 
These leaves Ac =^ F""^(c), c G M", are generically Lagrangean submanifolds of M. We shall 
always assume that F is proper, so that each Ac is compact. 

2.2 Known results 

Local and semi-global descriptions of completely integrable systems are provided by the following 
facts : 

• when the vector fields {Xi, . . . ,Xn) are independent at a point m, the Darboux-Caratheo- 
dory theorem states that the functions fi can be taken as "momentum" coordinates of a 
symplectic chart {(-x,^)} on a neighborhood U of m. In these coordinates, the foliation 
AcCiU for c close to F{m) is given by ^ = c. The Xi variables are therefore local coordinates 
for Ac. 

• when c is a regular value of -F - this includes the previous case, of course - the Arnold- 
Liouville theorem endows an invariant neighborhood Q, of Ac with action-angle coordi- 
nates : Q, is symplectomorphic to a neighborhood of the zero section of the cotangent 
bundle T*(T") of the n-torus = M^/Z", and the fibration F gets transformed into a 
smooth function of the momentum variables only. As a consequence, Ac, as well as the 
neighboring fibers, are Lagrangean tori equipped with an affine structure : the one given 
by the action which, in these coordinates, is linear. 

In the first case, the action is locally free at m, and gives rise to a local diffeomor- 
phism between (M'^,0) and (Ac, m). By transporting the standard Lebesgue measure via this 
diffeomorphism we obtain an invariant measure on a neighborhood of m in Ac which in Darboux- 
Caratheodory coordinates (x,^) is nothing else than \dx\. 

In the second case, this can be done globally to give a natural invariant measure pc on the 
whole torus Ac. Letting rric be the total mass of Ac, we can identify m~^pc with the Lebesgue 
(or Haar) measure on the torus T" via action-angle coordinates, pc smoothly depends on c and 
is called the Liouville measure. 
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Remark 2.1. This terminology sometimes also applies to the symplectic measure u)^^ /n\. 
Note that the push-forward of this measure by F gives a measure /ic on that satisfies 

= Pc® He, 

n\ 

and that is long known to be of particular importance when it comes to counting eigenvalues (see 
section |^). Because of the article [ pO[ [i^ is sometimes called Duistermaat-Heckman's measure. 

A 

Since we will be dealing only with pseudo-differential quantization, we will always assume 
that M is an open subset of a cotangent bundle T*X, which implies that the symplectic form 
is exact : uj = dao, where ao is the standard Liouville 1-form. If we let ic ■ Ac ^ M he the 
inclusion, the fact that Ac is Lagrangean implies that the 1-form i*ao is closed and thus gives 
rise to an element of the cohomology ff^(Ac,M). With slight abuse of notation, we will still call 
it [ao]- 

3 Semi-classical integrable systems 

The aim of this section is to define what we consider to be the proper semi-classical quantiza- 
tion of completely integrable systems (in accordance with ||]), and to present a geometric 
interpretation of their sub-principal terms. 

3.1 Definition 

Let X be a differentiable manifold of dimension n, equipped with a half-density Idxl^/^, and let 
Q be an open subset of T*X. The spaces of pseudo-differential operator that are used here are 
defined in section ^ however, we shall here only use the fact that they are defined up to the 
order 0{h?') by two functions of Q. : their principal and sub-principal symbols. 

Throughout this work, pseudo-differential operators are always classical and formally self- 
adjoint, so that their principal and sub-principal symbols are reaZ-valued. 

A set {Pi(/i), . . . , Pn{h)} of n pseudo-differential operators in \I''^(ri) of order zero is called a 
semi- classical integrable system if : 

• the principal symbols pi, . . . ,pn form a completely integrable system. 
. Vi,j [Pi{h),Pj{h)]=0. 

Note that the second condition already implies that the principal symbols are in involutions. 

Because of this definition any semi-classical integrable system has an underlying classical 
completely integrable system, which means that the main geometric ingredient of such a system 
is the Lagrangean fibration Ac given by the principal symbols. Of particular importance for us 
will be the principal Lagrangean Aq. 

Remark 3.1. All the results of sections 2-6 would still hold if the commutation property is 
weakened to [Pi{h), Pj{h)] = 0{h°°). Section 7 however requires the exact commutation of the 
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operators, in order to define their joint spectrum. Similarly, a number of results before that 
section would still be valid if the self-adjointness requirement is dropped and replaced by the 
assumption that both the principal and sub-principal symbols are real-valued. A 

3.2 The sub-principal form 

3.2.1 Definition 

Apart from the principal symbols, the data of a semi-classical integrable system defines another 
set of functions, namely the sub-principal symbols ri, . . . ,r„. We assume here that these func- 
tions are real-valued. While the principal symbols gave rise to the momentum map and its 
Lagrangean fibration, the sub-principal symbols will be viewed as characteristic of an infinitesi- 
mal deformation of the Lagrangean leaves. 

Definition 3.1 The sub-principal form Kc of a semi- classical integrable system is the differen- 
tiable 1-form on Ac defined at non-singular points of F by : 

Kci^i) = -Ti. 

Recall that Xi denote the Hamiltonian vector field of pj. At such a non-singular point, {Xi , ■ ■ ■ , Xn) 
is a basis of the tangent space of Ac. 

This form also appears in |^]. 

The first property of Hc is the fact that it is closed. Indeed, the r^'s are not just any 
functions. The Weyl rule (section |^, formula (^) applied to the commutation property of the 
pseudo-differential operators Pi is, at the order /i^, equivalent to : 

{n,Pj] = {rj.Pi}- 

In a symplectic chart given by the Darboux-Caratheodory theorem, this reads : 

dri drj 
dxj dxi ' 

and Kc is the closed 1-form Kc = — X] ^idxi. □ 

3.2.2 Deformation of Lagrangean submanifolds 

Let C{M) be the set of all Lagrangean submanifolds of a 2n dimensional symplectic manifold M. 
By the Darboux-Weinstein theorem, we can identify a tubular neighborhood of a Lagrangean 
submanifold Aq with a neighborhood of the zero section of r*Ao. Any other Lagrangean sub- 
manifold within that neighborhood can then be identified with a closed 1-form on Ag (see 
This gives a "chart" for £(M), making it formally a differentiable infinite dimensional manifold, 
whose tangent space at Aq is naturally identified with the space of closed 1-forms on Ag. 

One can be more specific. The exponential map for vector fields gives a diffeomorphism 
between a tubular neighborhood of Ag and the normal bundle -rf^ — . The latter is identified 
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with T*Aq by means of the symplectic form : let a; G Ag; to any X G T^M is associated the 
cotangent vector 

uj{X) = ixuJ\T,Ao- 

Since Aq is Lagrangean, the kernel of ui is exactly T^Aq. 

In this way, we can easily describe the tangent space of £(M) : an infinitesimal variation 

T\ A'l 

of Aq is by definition a vector field transversal to Aq, that is, a section of . By the above 
isomorphism, it can be identified with a 1-form on Aq. That infinitesimal variation is performed 
within the space of Lagrangean submanifolds if and only if this 1-form is closed, that is, the 
deformation vector field is locally Hamiltonian (see also @]). 

Let us apply this description to the case where a path in C{M) is given as a 1-parameter 
family of level sets of n independent functions p\ in involution : 

A* = {p\ = 0.} 



Lemma 3.1 The infinitesimal variation 1-form k = ^A*|j=o ^-^ given by 

^ dt \t=o 

Proof. Let x G A'^. Complete (Pi,---,Pn) i^^fo a Darboux-Caratheodory chart, so that the 
tangent vectors form a basis of the normal bundle. Taylor's formula ~ Pi~^i~§t\t—o ^^^^^ 
that at t = 0, the deformation vector field for A* is given by 







To this transversal vector field corresponds the 1-form k = ix^- On the basis (^pO, . . . , X^o) of 
TAO, it is given by K(;fpo) = u{X, X^o) = dp^.X = -^i^^q. □ 

Remark 3.2. For each t, the functions , . . . define a completely integrable system in a 
neighborhood of A*. The variation 1-form can be computed not only for A*, but also for the 

neighboring leaves A* =^ ni{pl)~'^{c). Then = ^A* |.j^q is given on Ac as in lemma : 



dp^- 

Kc{X o) = — i . 

at \t=Q 

The fact that for all c, Kc is closed can be checked by differentiating the equality = at 

t = Q. A 



Remark 3.3. It is natural now to be willing to consider the level sets 



as a "linearization" of the family A^. Unfortunately, these level sets are in general not La- 
grangean. More precisely, the symplectic form restricted to them is of order t^. However, one 
can prove that there exists a set of functions fi , . . . , f„ such that the functions 

~f def , j.~ 

Pi = Pi+ tri 

are indeed in involution, and such that the corresponding deformation 1-forms kc satisfies 

The value of fi on Ac is constant and obtained by averaging the sub-principal terms ^^^q 
over Ac- 

A 



3.2.3 Deformation of the Action integral 

Let A^ be a smooth family of Lagrangean submanifolds of M, and let ^t{9),0 € S*^ be a smooth 
family of loops such that each 74 is drawn on Aj. Suppose that on a neighborhood of the image 
of 7o, the symplectic form u is exact : lo = da. Then this holds for t small enough, and we can 
define the action integral : 

Mlt) = I a. 
Jit 



Lemma 3.2 The variation 1-form k = ^Atit=Q on Aq is characteristic of the infinitesimal 
variation of the action, in the following sense : 



d_ 
di 



Mlt)\t=o = K. 

J 10 



Proof. We want to prove that the 1-forms ^ {'yti^))\t=o ^^'^ 7o('^) cohomologous on S^. 
We have 



d_ 
di 



(7t(")) =7t* (^sjj^oA = 



= It [iBjj^da + diia-ao) ) . 
\ dt at / 

At t = 0, the vector field ^ splits into two components X^- and X^, the first one being tangent 

to Aq. The other one, in the normal bundle is by definition the deformation vector field 

of the family A^. Therefore, 

7o iiajtda) = 7o(^x,t^) + ll{ixM = + 7o(«)> 
\ at / 

which gives the result. □ 

Remark 3.4. Of course, the value t = Q plays a arbitrary role : if we define nt = ^Aj, then 
|A(7*) = K,. A 
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3.2.4 K as a semi-classical deformation 



Returning to the hypothesis of paragraph ( p.2.1 ), we are now able to give a geometrical inter- 
pretation of the sub-principal form. 

The total symbol Pi{h) of Pi{h) is considered as a semi-classical deformation of the principal 
symbol p^. Then around any regular point of the principal symbol, the sub-principal form is 
equal to the deformation 1-form - as h tends to zero - of any family of Lagrangean submanifolds 
defined as level sets of functions of the form + hvi + 0(/i^). 

Of course, the level sets of the total symbols themselves need no be Lagrangean. Recall 
however that in the case where is a regular value of the principal symbol we can average the 
sub-principal symbols without changing the cohomology class of Kc in such a way that the level 
sets of p^ + hri are indeed Lagrangean. 

As we shall see in section ^, the Bohr-Sommerfeld quantization conditions that apply to such 
a situation depend on the sub-principal symbols only through the cohomology class of Kc- This 
shows that articles like ^ that would rather assume that the sub-principal symbols are either 
equal to zero or at least constant on each Ac are not far from the general case. 

Another nice property of the cohomology class of Kc is its invariance under conjugation by 
unitary pseudo-differential operators. 

Proposition 3.3 Let Pi{h) he a semi-classical integrahle system with a non- singular principal 
Lagrangean Aq. Let V{h) he a classical unitary pseudo- differential operator on a neighhorhood of 
Aq. Let Pi = V~^PiV he the transformed integrahle system and k' the new sub-principal form. 

Then [n'] = [k]. 

Proof. Since P and P' have same principal symbols, we can write P' = P + hQ, with Q = 
{Qi, . . . , Qn)- The intertwining property VP' = PV reads : 

[P„V] = hVQj. 

At the principal level, this gives ^{pj,v} = vqj, which, writing v = e^'^, yields : 

{pj,c} = qj. 

The sub-principal symbols ^ of P*- ^ satisfy r'j = rj + qj, that is, 

K = K + dc, 

where dc denotes the differential of c as a function on Aq. □ 

Note that the unitarity of V is only used up to 0{h'^) in order to ensure that the sub-principal 
symbols of the transformed system are real-valued. 



4 Microlocal analysis 

The aim of this section is to present a minimal version of the microlocal analysis needed in order 
to give a precise and usable definition of what we call "microlocal solutions" of pseudo-differential 
systems (section |4.5| ) . 
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The results can be retrieved from [jl^, IC, 37, |3^; however, note that we focus here a little bit 
more on the microlocal point of view, and assumptions about "behavior at infinity" are usually 
irrelevant in this theory. 



4.1 Symbols 

We recall here the basic definitions of symbols and classical symbols. One reason for including 
here this standard material, besides the sake of completeness, is that our work mainly deals with 
local properties, in a neighborhood of a point or at least in a neighborhood of a compact region 
of the phase space; and for this purpose, a very simple notion of 'symbols' can be introduced, 
which is the following : 

Definition 4.1 A family of complex-valued functions {p{h))fi£H on a manifold Z , where H is a 
subset of M having as an accumulation point, is called a compactly supported symbol if there 
exists a compact set K C Z and an integer m £ Z such that every p{h) is a function with 
support in K such that : 

Va,3Ca, s.t. V/i G i/,sup|9°p(z;/i)| ^ C^/i'". (1) 

We denote by S^{Z) the vector space of all such symbols. If we drop the condition on the support 
we obtain the space of all symbols of order m, denoted by S"^{Z). 

The notation S^-^{Z) shall be used in any assertion that holds for SfP{Z) as well as for S'^{Z). 

Remark 4.1. In the notation S^^{Z), we forgot the dependence in H. Somewhat later we 
will have to be more careful, but this omission is harmless at this point. A 

Note that S'^^{Z) C SJq^{Z) for m ^ m' and, due to the Leibniz formula, the space 5*q-j(Z) = 
Umez '^{o)(^) symbols of any order is a graded algebra for the usual multiplication. It is also 
clear that SJ^^{Z) = If Z sits in a symplectic manifold, we can also endow 

with the Poisson bracket {piih),p2{h)} turning it into a graded Poisson algebra. 

Definition 4.2 A symbol {p{h))h^H £ ^^){^) said to be classical if it admits an asymptotic 
expansion of the form : 

in the sense that \/N ^ m, p{h) — J^k-mPk^'^ ^ '^(o)^^('^)' 

Remark 4.2. If we apply the above definitions to symbols independent of z, we get the notion 
of a "constant symbol". We will let Ch ^ S*{Z) denote the algebra of such constants. With 
that respect, the spaces S*{Z) and Sq{Z) can be considered as modules over C^. A 

A functional calculus can be performed on symbols of non-negative order; we will only need 
the following : if p{h) £ S'^^{Z) with m ^ 0, then expp(/i) G S^{Z). Conversely, let q{h) be 

an element of S9qJZ); if $7 is any proper simply-connected open subset of Z on which q{h) is 
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elliptic {i.e. \q{h)\ > a > on Q, for all h), then one can define a symbol p{h) G <Sq{Z) such 
that q{h) = ex.pp{h) on Q. 

Everything can be restated in the "classical" category. 

Definition 4.3 We say that a symbol is 0{h°°) if it is in S'^{Z) for all m > 0. We denote by 
S°°{Z) the subspace of all such symbols. Similarly, Sq°{Z) is the intersection of all Sq^{Z). 

Now, the spaces that will be really of interest for us are the quotient spaces S'^q,^{Z)/S^^{Z). 
Elements in the same class will be called microlocally equal on Z. This word is easily justified 
when Z is a subset of a T*X; its extension to any Z should raise no problem. The space 
of classical symbols modulo microlocal equality can then be isomorphically identified with the 
space of formal expansions of the form Y2k>TnPki^)^^ • Note also that the classical constant 
symbols modulo microlocal equality form the field of numerical formal series in h. In general, 
as it is easily verified, the quotient of Ch by microlocal equality is a field that we will denote by 

In order to get meaningful information from the "quantization conditions" that we are going 
to derive later, we need to consider symbols depending uniformly on a parameter E £ R". By 
uniformly, we mean that for every compact subset B of M", the estimate is valid uniformly 
for E G B. In the classical case, this is true for instance if each pk depends continuously on E. 

4.2 Pseudo-differential operators 

For the general theory, we refer to ITo). Proofs can usually be derived from the homogeneous 
theory of Hormander . We wish here to present a microlocal version of the standard classes of 
pseudo-differential operators. Let X be a differentiable manifold of dimension n, and M = T*X. 
In every local coordinates, a compactly supported pseudo-differential operator P{h) on an open 
subset U of X is an operator with smooth kernel of the form : 

^'^(^' y) = tAa^ I e^("-^)-«a(x, y, h)di, (2) 

where a{h) is a symbol in 5q([/ x f/ x M"). Note that has compact support in [/ x [/, hence 
these operators are continuous linear operators (for fixed h) from C°^(U) to Cg°(C/), or even 
from distributions to C^{U). We shall assume that a C°° half density \dx\'^ is given on 

X so that each function n on X is associated with the half density u|(ix|2. 

Pseudo-differential operators have a Weyl symbol, which is a /i-dependent function p{h) = 
aw{P{h)) on M, such that P{h) can be retrieved fromp(/i) by the so-called "Weyl quantization" 
scheme : 

{Op'^{p{h))u){x) = I e^(^-yy^p{^y,C;h)u{y)dydC (3) 

Now let p{h) G 5^(0), and x ^ ^^0° equal to 1 in some neighborhood of G M". We 
have a{x,y,^;h) = p{^^ , h)x{x — y) G S(P{U x U x M"), which allows us to form a pseudo- 
differential operator Op^ {p{h)) via (|2|). Then one can show that the Weyl symbol of the 
latter pseudo-differential operator, although perhaps not of compact support, is nevertheless 
microlocally equal to p{h) on Q. That motivates the following definitions : 
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Definition 4.4 Let he an open subset of M , in which we have symplectic coordinates {x,^,). 
We suppose here that x varies in U, an open subset of X containing tt{0,) where vr is the natural 
projection T*X ^ X . A compactly supported pseudo-differential operator P{h) on U will he 
called h-smoothing or "0{h°^)" on S7 or element 0/^*0° if its Weyl symbol belongs to the 
space 5°° (17). 

The quantization procedure Op^ {p{h)) is therefore independent of x modulo ^'o°(r2). 

Definition 4.5 We denote by ^"^{Q) the space of compactly supported pseudo-differential op- 
erators of the form OpY^ {p{h)) + ^'[^(17), for any p{h) and x defined as above. 



Remark 4.3. Dealing with pseudo-differential operators, we shall always assume that the 
symbols involved are classical. A 

This definition leads to an isomorphism between the spaces ^'™(r2)/*I'Q°(r2) and Sq^ {^l) / Sq° {0.) 
(for the latter is isomorphic to — ). This isomorphism, still denoted by aw, depends on the 
local coordinates, but we recall that the first two terms in the asymptotic expansions of the Weyl 
symbols are intrinsically defined on T*X, provided we let pseudo-differential operators act on 



half-densities (the proof of this was given in |21, prop 5.2.1] for homogeneous pseudo-differential 



operators, i.e. without a small parameter h, and applies to our situation with no essential 
change); these two terms are respectively called the principal and sub-principal symbols, and 
are compactly supported. 

It also ensures that theses classes of pseudo-differential operators are stable with respect to 
the operations of composition and taking adjoints. 

We are now able to define the notion of the restriction of a pseudo-differential operator; 
namely, if is any subset of fi, two elements P{h) and Q{h) of will be called microlocally 

equal on A^, and written "P ~ Q on iV" ifawiP) — ^^wiQ) on A^. Then we can perform inversion 
of elliptic operators in the following sense: P{h) G ^'J5^(J7) is elliptic at a point m if its principal 
symbol does not vanish at m. Then there is a neighborhood N oim and a Q{h) € ^'q '"(O) such 
that PQ ~ QP ~ / on N , where / is a pseudo-differential operator satisfying / ~ Op^(l) on 
N. 

Remark 4.4. The space of Weyl symbols is naturally equipped with the Lie algebra structure 
given by the symplectic Poisson bracket. On the other hand, the space of pseudo-differential 
operators has a natural Lie operator algebra structure. Though it is known that there is no hope 
for finding a "quantization" isomorphism that fully respects this Lie algebra structure, the Weyl 
quantization has the not-so-bad following behavior (see |2^) : Let pi{h) and P2{h) be symbols 
of order 0, and Pi{h) and P2{h) their Weyl-quantization. Then [Pi,i-*2] is a pseudo-differential 
operator of order 1 and of Weyl symbol 

-^[{pi,P2] + 0{h'')). (4) 

In particular the sub-principal symbol of [Pi, P2] is easily computed in terms of those of Pi and 
P2 ■■ 

0-sub(^[-pL,-P2]) = {Pl,0,O-sub(^2)} + {o-sub(-pL),P2,o}, 
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where pjfi is the principal symbol of Pj. A 

Remark 4.5. The above classes of operators would be purely abstract if we could not relate 
them to usual /i-admissible pseudo-differential operators in the sense of |37] and others, for which 
additional assumptions concerning the behavior of the symbols at infinity are needed. Such a 
link is here easy: if P{h) £ ^^(0,), where Q is some open subset of M", we can form a strongly 
admissible pseudo-differential operator on R" by extending the Weyl symbol to zero outside its 
support, and using a Op^ quantization. Conversely, if Q{h) is an admissible pseudo-differential 
operator in the sense of |3^, it has a Weyl symbol, that we can cut to make it compactly 
supported in $7, and get a P{h) again. If Q{h) was already obtained by a P{h) in ^'0^(17), then 
by definition the new P{h) will microlocally coincide with P{h) on the set where the cut-offs 
are both equal to 1. The interest of this microlocalization procedure will be made clear in the 
next section (prop [4. 3D . A 



Finally, compactly supported pseudo-differential operators on manifolds can then be defined 
to be locally finite sums of operators admitting the above description in a local coordinate chart. 
The principal and sub-principal symbols remain well-defined functions on T*X. 



4.3 Fourier integral operators 

Similarly to the space ^'q(O), microlocal classes of Fourier integral operators can be constructed. 
We will here skip most of the details. Let X, Y be differentiable manifolds, and ^Ix and Qy 
open subsets of T*X and T*Y respectively. Suppose we have a symplectic diffeomorphism x 
from fiy to Qx- Let A be the graph of x in x ^Y- It is an immersed Lagrangean submanifold 
of T*X X T*Y, where T*Y denotes the cotangent T*Y equipped with —lo, the opposite of the 
standard symplectic form. Thus, A admit parameterization by non-degenerate phase functions 
ip : A = UA^. For details about this construction, refer to [|l8[ or Q. Following |jlO|, we define 
^Q^{i}x,^Y,x) to be the space of operators whose kernel is microlocally equal to a compactly 
supported classical /i-oscillatory integral on A (see paragraph [4.61 ). On each A^, it has a principal 
symbol defined as a function on A^. A Fourier integral operator U{h) is said to be elliptic at 
a point {x,(,) or, equivalently, at the point x~^{^jOj if its principal symbol does not vanish at 

{x-Hx,0,ix,0) e K- 

These operators behave as expected with respect to composition and microlocal inversion of 
elliptic operators. Most important for us, they allow to find normal forms for pseudo-differential 
operators thanks to the celebrated theorem of Egorov : if P{h) is a pseudo-differential operator 
on ilx of principal symbol p and U (h) an Fourier integral operator elliptic at a point (x, then 
U~^PU is, near x~^(^)0) ^ pseudo-differential operator of principal symbol po x- 



4.4 /i-admissible functionals 

In the above, we have let pseudo-differential operators act on /i-independent distribution half- 
densities, which is of course not enough for our purposes. Let U he a regular domain in M". 
If u £ V'iUjQi), then for any pseudo-differential operator P(h), Pu is a compactly supported 

2 

distribution, hence in the Sobolev space for some integer s; this holds for any fixed h. In 
order to get asymptotic information for all h £ H, we are lead to the following definition (see 
13^, prop IV-8]) : 
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Definition 4.6 A family {uh)h,zH of distribution half densities is called admissible on Q if for 
any pseudo- differential operator P{h) G ^'g(ri), there is an N & 2, and some s G Z such that 
h^{Puh) is uniformly bounded in H~'^{U) for all h £ H. We denote by Vfi{Q) this space of 
admissible functionals. 

Dfi{0,) is then a C/i-module that is by definition stable under the action of ^1^(0,). 

Proposition 4.1 If Uh is admissible on Q then for any P{h) G ^o(^^) there is an N' ^ Z, such 
that h^ Pufi is bounded in L?'{U). 

Proof (see p. 195-196]). One can find an elliptic operator D{h), that is a zero-order 
usual pseudo-differential operator on U, that uniformly maps L'^{U) into some power of h times 
H~%U) (extend to M", and take D = /i^(-A + |xp + ly/'^). Let D'^ be a right inverse 
of D, and P{h) G ^o(^^) be as in definition So there is a N' and a C < co such that 
\\D-'^{h)P{h)uh\\L2{u) < Ch~^' . Now let il' D O and let I{h) G have symbol equal 

to 1 on a neighborhood of $7. Then (/(/i) — l)P(/i) G ^'g°(r2), and by the symbolic calculus, 
I{h)D{h) G Therefore, I{h)D{h) is uniformly L^-continuous th.II-36]), and we can 

write \\P{h)uh\\L2 = \\I{h)D{h)D-^{h)P{h)uh\\L^ + 0{h°°) ^ C'h^^' . □ 

A natural notion of "microlocal equality" for admissible functionals can now be defined : 

Definition 4.7 T)'^{VL) denotes the space of Uh G Vhi^) such that for any pseudo-differential 
operator P{h) G ^o(^) have : 

WPuhWmu) = o{h°-). 

Two admissible functionals and Vh are called microlocally equal on J7 if they belong to the 
same class modulo Following JI7[ /, we will write "uh ~ Vh on il" in that case. If 

m G T*X, we say that Uh ~ Vh at the point m if there exists an open neighborhood Q of m such 
that Uh ^ Vh on Q. 

Note that, because the above definition only involves estimates on compact subsets of il, Uh ~ Vh 
on 0, if and only if Uh ~ Vh at each point of 0,. Moreover, to test microlocal equality at a point, 
it is sufficient to pick up an elliptic operator that satisfies the required estimate. There again 
we have obvious although useful properties : 

Proposition 4.2 • if P{h) G ^"^{Q), then for all Uh G Vh{'^), Put ~ on 0; 
• if Uh ^ on $7, then for all P{h) G ^'o(^)' ^'^h ~ on $7. 

This is essentially due to the L^-continuity of pseudo-differential operator of order 0. In other 
words, the quotient space ^q{VL)/'^'^{^) acts naturally on P/j(ri)/P^(il). This says that the 
action of a pseudo-differential operator on an admissible functional is microlocally given by 
the (formal) Weyl symbol ayi/. Note also that the space Vh{^) /V^{VL) has the structure of a 
C/j-vector space. 

In particular, if is any subset of $7, the definition gives a natural microlocal equality on 
A, which is compatible with the notion of restriction of a pseudo-differential operator previously 
introduced: two pseudo-differential operators microlocally equal on A have the same action on 
two admissible functionals microlocally equal on A. 

The same remark applies for the following interesting result : 
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Proposition 4.3 If Q{h) is a global admissible pseudo-differential operator in M" in the sense 
of W\] j ^ o-'n- open subset of M^", Uh an admissible functional in O, and P{h) S ^o(f^) a 
microlocalization of Q in a compact K ofQ, then Quh is admissible in Q and : 

Quh ~ Puh on K. 

Alternatively, one can relate this notion of microlocal equality for admissible functionals to 
the so-called semi- classical wave front set, as introduced in |^], and denoted by WFh{uh)- 
This notion is not stabilized yet, in the sense that, depending on authors, it includes or not 



uniform estimates at infinity (see remark in |11, p. 1541]). Anyhow, the simpler part of WFh 



that is to say its intersection with T*X, is defined as follows : 

suppose Uh is admissible on T*X. Then WFh{uh) n T*X is the complement subset in T*X 
of the biggest open subset Q such that ~ on 17. 

The definition is of course based upon the homogeneous wave front set introduced by 



Hormander (see |3C, sec.2.5]). It has a useful local characterization in terms of the semi-classical 
Fourier transform : 

Lemma 4.4 (|37|, [1C)[ |) An admissible functional microlocally vanishes at the point {xo,^o) £ 
T*X (equivalently, (xo,Co) WFh{uh)) if and only if there exists a function ip G C^{X) with 
ip{xo) 7^ such that, in some local coordinates : 

:Fhiipuhm = 0{h°^) 

uniformly for in a neighborhood of ^o- 

Here J-'h denotes the "semi-classical Fourier transform" : 

(27r/l) 2 JiRn 

(we will sometimes use the usual Fourier transform T = T\). Note that Th is a priori defined 
for tempered distributions, but the Fourier transform of a compactly supported distribution in 
can be thought as a continuous function (it is analytic). For the proof of this lemma, one 
constructs a compactly supported pseudo-differential operator with "rectangular support" , that 
is with symbol of the form x(Ov'(^)- 

Unfortunately, because of the lack of uniform estimates in the fibers one should not expect 
an analogue of theorem 2.5.3 of the previously cited article [^]. Indeed, if U is an open set of 
X, then WFh(uh) n T*\J = does not imply that is locally 0{h°°) on U. For instance, one 



easily sees from lemma iA that Uh{x) = a{x)e^^^^ is microlocally zero on T*M", while obviously 
not locally 0(/i° 



The solution found in |10| for functionals on is to extend the wave front set to a subset 
of the compactified cotangent bundle obtained by adding a point at infinity to each direction ^. 
An admissible functional Uh is then said to be microlocally zero in the direction (2;o,M'*'^o) if; 
for large ^ in a conic neighborhood of and for any cut-off function if, the Fourier transform : 

is of order 0(/i^/|^|^) for all > 0. Then the following statement is valid : 
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Lemma 4.5 An admissible functional Uh is 0{h°°) uniformly for x near Xq if and only if there 
exists a function ip G Cq°(X) with ^{xq) 7^ such that : 

for all(,£ W\ for all N £N. 



The condition expressed by this lemma is that there exists a neighborhood U of Xo such that Uh 
is microlocally zero in T*U as well as in all directions of T*U. 



Similarly to lemma ^4.4| , we have for instance the following fact : let Uh be an admissible 
functional defined on M", so that J-hUh makes sense. If J-hUh is a 0{h°°) function in a neigh- 
borhood of any ^ in an open set T, then is microlocally zero in the whole x T. This 
allows the following construction : let 1 = xi + ' ' ' + Xfc be a partition of unity subordinated to 
the open cover |J^. Tj in the ^-space. Let U be an open subset of M" with compact closure. If 
Vj,h, j = 1, ■ ■ ■ ,k are admissible functionals on ?7 x Tj that microlocally coincide on non-empty 
intersections U x Tj^ n • • • n Tj^ , one can define an admissible functional on U x \Jj Tj by 
cutting off the fj,h's outside U in M" and requiring : 

k 

Then Uh microlocally coincides with each Vj^h on U x Tj. Indeed, from the microlocal equality 
of the Vj^s on non-empty intersections, and the compactness of their support, it follows from 
lemma [4.4| that the Fourier transforms J-hVj.h are C°° functions that coincide with each other up 
to 0{h°°) around each point of such non-empty intersections. Thus J-hUh also has that common 
value around these points, which yields the result. □ 

Finally, if U{h) G ^'q'{VIx-,^YiX)i then for any admissible Uh in Qy, Uuh is admissible in 
Vtxi and because of Egorov's theorem, U transports the wave front set by x- 



4.5 Microlocal solutions 

Definition 4.8 Let i7 he an open set in T*X . An admissible functional G Dhi^) is called a 
microlocal solution on of a system (Pi(/i), . . . , Pk{h)) of pseudo-differential operators if : 

yj = 1, . . . ,k, PjUfi ~ on 

We allow here to choose a suitable set H where h is to vary. The space of all microlocal solutions 
of {Pi, . . . , Pk) modulo microlocal equality will be viewed as a C/i-vector space. Because of 
microlocal inversion of elliptic operators, any such solution is "microlocalized" in (~]jPj^{0), 
which means that WFfi{ufi) must lie in this set {pj denotes the principal symbol of Pj). If two 
solutions have different non-empty wave front sets, then they are independent. In case the pj's 
are real- valued, the "propagation of singularities" theorem states that WFh{uh) is locally stable 
under the action of the Hamiltonian flow of pj. 

As we will see later, global problems often have no solution in the previous sense. Suppose 
that the Pj depend continuously on an additional parameter in a compact topological space 
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V (cf. end of section 4J). Depending on the degeneracy of the p/s, it might be impossible to 
find a solution (in the previous sense) for some E, and depending on E, the appropriate sets He 
where h varies might be disjoint. Still, the whole theory being asymptotic, we need to control 
things when h tends to zero ! Therefore, we introduce the following : 

Definition 4.9 Let T be a subset o/]0, 1] x V. Let H be the projection of V onto ]0, 1]. We 
assume that is an accumulation point of H . We say that a family {u(^h^E)){h,E)eT ^■^ admissible 



if definition |^.q holds uniformly for {h, E) G T . 



Recall that E varies in a compact. Accordingly, 'U{/i,_e) microlocally vanishes at point m if there is 
an elliptic pseudo-differential operator P{h) at m such that \\Pu(^ii^e)\\l'^{u) — 0{h°°) uniformly 
for {h,E) G r. Thus we are still able to define a semi-classical wave-front set WFh{u^h^E))- 

Now if U(j^^E) is admissible and P^" depend continuously on E, then P^Uf^j^ ^) is admissible, 
and we say that U(/i,_e) is a microlocal solution of the system {Pf , . . . , P^) on Q if P^U(^fi E) ~ 
on 17 according to that definition. Then we have : 

Proposition 4.6 Let u^^i^^e) be a microlocal solution of the equation P^u^^f^ E) ~ on (7. The 
following inclusion holds : 

WFh{u(^f,^E)) C {{x,0 G T*X, 3Eo G V^, PeAx,0 = 0}, 

where Vq is the set of Eq G V such that {0,Eo) is an accumulation point ofV. pE is the principal 
symbol of P^{h). 

Proof. Let {x,(,) G T*X such that for any Eq G Vq, the principal symbol pEo is non zero at 
(x,^). Since Vq is compact and pE is continuous in E, then there is an open neighborhood V of 
Vq and a constant c > such that for E £ V, \pe\ > c. Then there is an open neighborhood Q of 
(x, ^) on which that remains valid. The complement being a compact set of non-accumulation 
points, we can pick up a ho such that ([0, ho] x V'^) H F is empty. From now on, we will therefore 
restrict F to [0, ho] x V. Because of the uniform ellipticity of P^ on Q, for (h, E) in that new F, 
we can construct a pseudo-differential operator Q^{h) depending continuously on £' G V such 
that Q^P^ ~ / on $1. It is not difficult to see that if an admissible V(^fi,E) is microlocally zero 
at a point m, then for any pseudo-differential operator Q^{h) depending continuously on E in 
some compact, then Q^v^^^e) is still microlocally zero. Back in our problem, we see that : 

U{h,E) ~ at (x,^), 

which is what we needed to prove. □ 



4.6 Oscillatory integrals 

The WKB ansatz for solving Schrodinger-type equations consists of restricting the scope of 
solutions to a subspace of admissible functionals : the space of oscillatory integrals, that we 
wish to describe now. These functions have a long history, particularly in Quantum mechanics, 
and it is well known since the treatise of Maslov |^^, made more rigorous and expounded by 
Hormander and Duistermaat (see in particular ||l8| ) , that they are particularly fit to locally solve 
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generic microlocal pseudo-differential equations. In fact, as we sliall see later, in some cases all 
solutions are microlocally equal to oscillatory integrals. 

We recall here the definitions of []T8| , with some precisions. 

We are given an compact immersed Lagrangean manifold A G T*X, endowed with a smooth 
half-density p. Suppose that is fixed a covering of A with simply connected embedded open 
subsets ^.^p^ described by reduced phase functions ipk{x,9) ("reduced" means that the number 
of additional oscillatory variables 9 is exactly the maximum dimension of the kernel of dn, 
where vr is the projection A^, X. This dimension is known to be the minimum dimension of 



variables needed for a phase function defining A - [30, theo 3.1.4] or [18, 1.3.6]). (f gives rise to 
a function on A^ - that we still denote hy ip - whose exterior differential is the Liouville 1-form 

OiO- 

An oscillatory symbol on A<^ is a half-density cr(/i) on A^p of the form : 

a{X;h) = c{h)ei'^^^^a{X;h)p{X), 

with c{h) G C/i, and a{h) a classical symbol in S*{Aip). Will will say that a is classical if c is. 

If a{h) is of order and elliptic on A^ (as we will generally assume), we will prefer the 
notation : 

ct(A;/i) = e*'=('^)+^*(^;'^V(A), (5) 

where c{h) G Ch, and $(/i) is a iS^-valued classical symbol on of order —1, whose principal 
symbol satisfies : 



As usual, we denote by the local diffeomorphism from Cip to A,^ given by {x, 9) i— > 
(x, If (x,e)), where = {{x,9) eXx rN,^{x,9) = 0}. 

An oscillatory integral on A^ is a half-density on X of the form : 

Uh{x) = ^^^^ ^ I eT^^^''''^h{x,e;h)d9\dx\^, (6) 
{2i-Kh)- Jr^ 

where c(/i) G C/^, and b{h) is a classical symbol in Sq{X x M^). Then 

is an oscillating symbol. Uh is an admissible functional, and because of the stationary phase 
formula, is uniquely defined modulo microlocal equality by the asymptotic expansion of b{h) on 
any neighborhood of C^. a is thus well defined by Uh modulo microlocal equality. 

Definition 4.10 An admissible functional is a Lagrangean distribution on A if there is a 
cover A = UA^ such that Uh is microlocally equal to an oscillatory integral on each A^. 

On each A^, such Lagrangean distributions have a symbol a{h) defined modulo 0{h°°), but in 
general {N ^ 0) the morphism ^ c(/i) is not injective; however, one can prove that it is 
an isomorphism at the principal level (if Uh can be written with b{h) of order m such that the 
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h"^ part of (t(/i) vanishes on a neighborhood of {xo,Co) £ A^, then near j^^(xo,Co)) Uh can be 
rewritten with a of order m + 1). 

On the other way round, the symbol map is surjective, and one can construct right inverses 
("Maslov canonical operators") in the following way : let o"(/i) be a oscillatory symbol near 
ixo,(,o) £ A^. Then it defines the function b{h) on a neighborhood of bo = j^^{xo,(,o) in C^. 
If iV > 0, is a proper submanifold of X x M^. Let Bo be a sub-bundle of T{X x R^) 
supplementary to TC^ near bo- One defines the germs of b on C^p to be zero along Bq^ and 
smoothly extends 6 to X x accordingly. Then b has the same /i-order as cj(/i) and the 
oscillatory integral Uh = Op{a) so constructed has symbol cr{h) (mod 0{h°")). 

Of course, since we have not made any mention to the Keller-Maslov bundle yet, the symbol of 
a Lagrangean distribution depends on the phase function ip considered. The result of Hormander- 
Maslov is the following : if and Vh are zeroth order oscillatory integrals on and 
respectively that are equal modulo 0{h) at a point (x,^) G n A^, then there is an integer jj, 
such that, at the principal level, cJu^ ~ e^^^a^i^ + 0{h) at (x,^). Because of the normalization 
in @, // = //^ — //<^, where, for a phase function g, 2^g is the signature of the 0-Hessian of g at 
(x, ^) minus Ng. 

With these transition functions, one can globally define the principal symbol of a Lagrangean 
distribution as a section of the Keller-Maslov bundle. In the sequel, Maslov indices will always 
be prompted very clearly, so we don't insist any further on that point. 

Note that if Uh serves to construct a Fourier integral operator U{h) as mentioned in section 



4.3| , then the principal symbol of Uh is invariantly defined by U. 

If the symbol is elliptic and we write it as @, we will refer to ^{h) as the total phase of Uh, 
and the principal phase is -|- <I>o, both defined modulo a constant c{h) G C^. 

Finally, we admit the dependence of a Lagrangean distribution on an additional parameter 
E, provided the involved estimates are uniform in the sense of definition 4^ . It is the case for 
instance if the symbol b depends continuously on E. 



5 Non-singular quantization conditions 

Let Pi{h), . . . ,Pn{h) be a semi-classical integrable system, and denote by Aq the principal La- 
grangean manifold associated to it. Aq may be critical, but in this section we shall work on a sub- 
manifold A C Aq consisting only of non-singular points of the momentum mapping (pi, . . . 
Moreover, we shall always assume that A is compact. 

With these data, we express a necessary and sufficient condition for the existence of mi- 
crolocal solutions of the system on A (proposition ^I^ ). Thanks to the WKB construction, we 
prove that this condition involves, at least at the principal level, geometrical characteristics of 



the system, namely the Liouville 1-form ao and the sub-principal form k (theorem 5J 

As an application, we recover the well-known Maslov-Bohr-Sommerfeld quantization condi- 
tions for 0{h^) solutions in the case where A = Aq is a whole non-singular Liouville torus. 

When the system depends smoothly on an n-dimensional parameter ii^ in a suitable sense, 
the above results are seen to be uniform on E, which will allow us to apply the quantization 
condition to the study of semi-classical spectra. In the case of a non-singular Liouville torus, 
the study of the spectrum was also completely carried out, with a different approach, in ||5|. 
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5.1 Dimension of microlocal solutions 

The following result is a generalization of lemma 18 of [0]. 

Proposition 5.1 Let x £ A be a non- singular point o/Aq. The space of microlocal solutions of 
the system PiU ^ Q on a neighborhood of x, modulo microlocal equality, has dimension 1. 

If X Aq, the space of solutions of this system has dimension 0. 
Proof. The second point comes from the fact that near x, the system is elliptic, hence microlo- 



cally invertible (see proposition 4.6). 



Concerning the first one, we use the following microlocal normal form, which is standard : 

Lemma 5.2 There exists an elliptic Fourier integral operator U{h) associated to a Darboux- 
Caratheodory coordinate chart (x,^) such that 

I OXj 

on a neighborhood of x. If every Pj is formally self-adjoint, then U can be chosen to be unitary 
(mod 0{h°°)). 

Accordingly transformed, the system always admits the constant solutions. It remains now to 
prove that these are the only ones, that is, that any solution Uh of -£pu ~ is microlocally 
constant. Let Uh be such a solution. We can always assume that it is compactly supported. We 
want now to "cut-off in the frequencies" , which we do in the following way : 

Let Vh an admissible functional defined by J^h^h = X^hUh, where the cut-off function x 
is and has value 1 on a neighborhood of 0. Then Uh and Vh are microlocally equal on a 
neighborhood of 0. Moreover, is a local solution (modulo 0(/i°°)) of the system near 0. This 
comes from the fact that ^h{j^^) = ij^hVh satisfies the hypothesis of lemma 4.5 . 

We have thus reduced the problem to a system of the form : 

dhVh = wh = 0{h^) 

on a neighborhood of (here dh denotes the exterior differential operator d^ = (fgf^, • • • , fgf"))- 
The solutions of this system are necessarily of the form 



CStCh + T Wh(tx).xdt. 
h Jo 



The second term in the sum is 0{h^), as Wh was. Hence it is microlocally zero, which gives the 
result. □ 

Let ipt = V'i(ti) o • • • o (/7„(t„) be the classical flow associated to the momentum map of the 
system. By transport of singularities, we obtain the following corollary : 

Corollary 5.3 // the action of (pt on A is transitive, then the space of microlocal solutions of 
the system PiU ~ on A has dimension at most 1. 
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If we restrict to a simply connected open subset of Aq, and if all trajectories of the flow 
remain open (which is automatic if fl contains no singular point), then local microlocal solutions 
can be glued together to obtain a non-trivial solution of the system on J7. This occurs for instance 
if 0, is obtained by propagating a non-singular point of Aq by the flow during a short time (cf. 
theorem 1.4.1 of 

As soon as we have closed trajectories, we get obstructions, which we describe in the following 
paragraph. Later on (paragraph ( ^.4D ), we will show how to interpret these obstructions thanks 
to the WKB method, and get the so-called "Bohr-Sommerfeld quantization conditions". 

5.2 The bundle of microlocal solutions 

For X G A, let Ch{x) be the 1-dimensional vector space over Ch of microlocal solutions of the 
system at the point x (modulo microlocal equality). For any open subset Vt of A, the disjoint 
union Ch{^) of the Ch{x), x G Q has a natural structure of a C/i-line bundle over Q, defined 
by the following local trivializations : if = Ufio,, where every is chosen small enough so 
that the system admits a non-trivial microlocal solution on it, then we define (pQ^ to be the 
isomorphism : 



{x,c{h))^c{h)ul (7) 

If we restrict to non-trivial solutions, this also defines a principal bundle >C/i(A) over il. with 
structure group C^* = \ {0}. This bundle has a natural flat connection, whose parallel 
sections are those that are locally constant in the above trivializations (from the infinitesimal 
viewpoint, its horizontal distribution is defined to be the image of the standard horizontal spaces 
of X C/i by the trivialization maps v^n^)- It is well defined because on intersections H i7/3, 



such trivializations differ by a constant, due to proposition 5.1. By the same argument, it is 
also independent of the set of microlocal solutions chosen to define the trivializations. The 
interest of this connection is stressed by the following lemma : 

Lemma 5.4 There is a natural identification of the space of non-trivial microlocal solutions of 
the system on an open subset Q C A with parallel sections of Ch{^)- 

Proof. Let us first consider the case = ri^. The general case will follow by a microlocal 
partition of unity. Now, using the trivialization Lp^a-, a section of Ch{VLa) takes the form : 

ah{x) = c{x]h)ul. 

This section is parallel if and only if c(x; h) = c{h) is a constant (modulo microlocal equality. 



of course), so, by proposition 5.1, if and only if (^^{x) = Gh'is & microlocal solution on Q^- It is 
non-trivial if and only if Ch 7^ 0. 

Now assume is any open subset of A. One sense is obvious : if Uh is a non-trivial solution 
on 0, then it can be used to trivialize the line bundle Ch{Q) as in (^). Now, since the definition 
of the connection is independent of the chosen solutions, the constant section ah{x) = Uh is a 
parallel section. 
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On the other way round, suppose crh{x) is a parahel section of Cfi{Q). Then we proved 
that its restrictions to each 0,^ are non-trivial microlocal solutions on Oq. Let P'^{h) be a 
pseudo-differential partition of unity on O subordinated to the cover 0.^ (it rneans that the Weyl 
symbols form a Cg° partition of unity of Q). Then each P'^v^ is an admissible functional on a 
neighborhood of fi, and the sum 

a 

is microlocally equal to (and hence to af^) on each ^1^- Therefore, Vh is a microlocal solution 
on fi, whose associated parallel section is af^. This terminates the proof of the lemma. □ 

Thus, the obstruction to the existence of a microlocal solution on A only lies in the holonomy 
of the connection, which, in this case, describes the monodromy of the solution. This holonomy 
can be represented by the cohomology class of a Cech 1-cocycle on A with values in C/j*. 
Summarizing, we get : 

Proposition 5.5 There exists a non-trivial microlocal solution of the system on A if and only 
if the holonomy cocycle of the -principal bundle Ch{-^) with its natural fiat connection has 
microlocally trivial cohomology. 

From now on, our principal task will be to explicit that holonomy. Besides the above definition, 
we can give a direct and more concrete description of it in terms of Cech cohomology, that shall 
be more efficient for our further analysis. 

We still let Oq, be a (finite) cover of A by simply connected open subsets such that each 
0,a admits a neighborhood in M on which proposition applies. We thus get a family tt^ 
of non-trivial microlocal solutions of the system on fio,. Moreover, on non-empty intersections 
Qa n 0^, the same proposition affirms the existence of unique non-zero constants c^°'{h) G 
such that : 

~ on n fif^. 

The c^"'s satisfy the cocycle condition : 

on every non-empty intersection Oq, H n 17^. Furthermore, another choice of solutions u"f^ 
would give rise to constants c'^° differing by a coboundary d°'{h) : 

ifSa d^ 

d°' 

This Cech cocycle c^" is exactly the holonomy form of the connection introduced above. The 
matching conditions are then the following : 

"there exists a non-trivial microlocal solution on A if and only if c^"(/i) is a coboundary 
(modulo 0(/i°°))." 

It is, of course, a simple re-formulation of proposition In the sequel, instead of the 
complex- valued Cech cocycle [c(/i)], we will rather deal with the two real- valued cocycles obtained 
via the exponential representation : 
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in which p^'^{h) € M and A^°(/i) G T. We denote by p{h) and X{h) the corresponding Cech 
cocycles. 

Proposition 5.6 Suppose that every Pj is formally self-adjoint. Then [p{h)] is microlocally 
trivial, and [A(/i)] admits an asymptotic expansion of the form : 

Proof . We cover A by open subsets Via on which the normal form of lemma p.2| applies. 
This gives rise to a family of microlocally unitary Fourier integral operators C/"(/i) on Oq,. An 
associated set of solutions is given by = C/"la, where Iq, denotes an admissible functional 
microlocally equal to the constant 1 on 17^. Suppose JIq H 17^ is a non-empty intersection and 
let X be an element of it. Finally, let Ix denote an admissible functional microlocally equal to 
1 on a neighborhood of x. On x, we have uf^ ~ {7°lx and ~ U^lx- Therefore, the constant 
c^'^{h) is given by : 

The unitarity of the Fourier integral operators implies that |c^"| ~ 1, which gives the first 
statement of the proposition. 

Now, (C/°)-^f/^ is an elliptic Fourier integral operator which implies that (C/") ^U^lx is 
at X a classical Lagrangean distribution with a symbol of the form (^). This gives the last 
statement of the proposition. □ 

Remark 5.1. The description of the cocycle c^"(/i) involves no regularity on A. And indeed, 
the notion of a flat line bundle with holonomy is well-defined on any topological manifold. This 
will enable us, in section |6.7| , to apply that description on a singular Lagrangean. For the 
moment however, the C°° structure of A makes it possible to apply the De Rham isomorphism 
and get a smooth connection 1-form \{h) on A. In the next sections, precisely, we show how the 
1-forms ao and k enter the picture. A 



5.3 WKB method 



Although this might not seem obvious, propositions 5.5 and 5.6 are essentially related to the 
celebrated WKB construction. 

The WKB method for solving pseudo-differential equations consists in looking for a special 
kind of approximate solutions, called oscillatory integrals (see section The use of semi- 

classical oscillatory integral in mathematics was initiated by Maslov, and pursued by different 
authors including Arnold, Hormander and Duistermaat. The result is that, provided we are 
given an invariant half-density on A (this will come from the Liouville measure), we can always 
construct such solutions in a neighborhood of a non-singular point of Aq. According to propo- 



sition 5.1, such a solution has to be the only one - up to multiplicative constant. This means 
that we can freely use the WKB construction to look for general properties of our microlocal 
integrable systems. This is the goal of the following paragraph. 
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Note that we previously constructed microlocal solutions via another method, namely the 
use of a normal form. The link with oscillatory integrals lies in the use of a Fourier integral 



operator; more precisely, with the notations of lemma 5.2, the admissible functional Ul is a 
Lagrangean distribution that is a microlocal solution of the system on a neighborhood of x (1 
denotes a functional microlocally equal to the constant 1 on a neighborhood of x). 

The derivation of the transport equations in the WKB method for completely integrable 
systems is fairly standard now (see for instance [^), and our aim here is only to briefly recall 
the idea. Of particular importance is the entrance of the sub-principal form n into the picture. 

Let Uh be a Lagrangean distribution on Aq of order 0. To ensure that it is not trivial, we 
assume that it is elliptic at a point (xq, ^o) S Aq, so that we use the notation (||) for the symbol 
cr{h) (we assume that we work on a fixed A^). The viewpoint here is to find conditions on a in 
order that Uh be a microlocal solution of the completely integrable system at (xo,Co)- 

First of all, if one looks for a solution up to 0{h) terms, one finds the condition that is 
indeed a phase function for A<^, i.e. : on A^, 

h " T' 

Note that this first condition only involves the principal symbols of the Pj^s. 

Keeping <I>_i fixed, one can look now for a solution up to 0{h?). Since we suppose that the 
sub-principal symbols are real (which occurs for instance if the Pj's are formally self-adjoint), 
then this raises two necessary and sufficient conditions : a) the half-density p must be invariant 
with respect to the Hamiltonian action of the system; b) on A<^, one must have 

where n is the sub-principal form of the system. Because of condition a), we can assume that p is 
the canonical Liouville half-density on Aq. Condition b) implies that k must be closed, which is 
already given by the fact that the Pj's commute with each other up to the order 0{h?). Finally, 
note that if the sub-principal symbols are complex- valued, the conditions are still sufficient, but 
not necessary. 

The previous step generalizes by induction to provide solutions up to any order. However, 
notice that the two conditions obtained above depend only on the principal symbol of Uh and 
on some closed 1-forms uniquely defined by the system. The next conditions involving higher 
order terms of the symbol won't of course be invariantly defined anymore. Still, there is a way 
to keep track of the total symbol of u^, and this is to fix a "Canonical operator" (see section 
Uh = Op{a{h)), as stated in the following proposition : 

Proposition 5.7 Let he a neighborhood of (xo,Co) described by a phase function 

if. Let Op be a Canonical operator on A,^. There exists closed real-valued 1-forms k/., k = 
— 1,0, 1, . . . , oo on Aq such that : 

The oscillatory integral Uh = Op{e^^^^''^'' p{\)) is a microlocal solution of the system up to 
the order 0{h^) on A^, if and only if the ^k's, k = —1, . . . , N — 2 satisfy : 

d^k = Kk- 

At_i is the Liouville form oq restricted to Aq, and kq is the sub-principal form k. For k > 0, 
Kfc depends on the Op-quantization. 
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A similar result holds if the Pj's depend continuously on a parameter E. Note that in ||8|, 
p. 38] Y.Colin de Verdiere included the dependence on a parameter directly in the definition of 
the oscillatory symbols. 



5.4 The regular quantization conditions 

Theorem 5.8 We still assume that the Pj 's are formally self-adjoint, and that A is a regular 
open subset of Aq. 

Let [X{h)] G H^{A,T) be, as in proposition |5.4 the holonomy phase of the microlocal bundle 
Cfi{A). There exists a non-trivial microlocal solution of the system on A if and only if 

[X{h)] G H\A,2TTZ) + 0{h'^), 

and any other solution is then a multiple of it. Moreover, [A] admits the following expansion : 

oo 

k=-i 

where fi is the Maslov cocycle, A_i is the restriction of Liouville 1-form uq on A, and Aq = k^a 
is the sub-principal form of the system. 

Proof. The first points are direct consequences of propositions 5^, 5^ and |5.1| . It remains 
to prove the statement about the principal terms of the asymptotic expansion of [A]. We cover 
A by simply connected open subsets Oq, such that each Oq, admits a representation by a non 
degenerate phase function ipa. The most immediate consequence of proposition 5.7 is that 
microlocal solutions as Lagrangean distributions always exist on Q^; Let us denote by such a 
set of solutions, and by their principal phases. Finally, we let fi"^ be the (integral) value 

of Maslov's cocycle associated with that cover. 

The equation 

u1 ~ c^''{h)ui on 

implies that argc^" is, mod 0{h), the difference of the principal phases of and u^, plus the 
value of the corresponding Maslov cocycle. That is, 

argc^"(/i) = - + //"/^^ + 0{h). 

Because of proposition 5.7, $° — is the value of a Cech cocycle whose cohomology is repre- 
sented, via the de Rham isomorphism, by the cohomology of the smooth 1-form + k restricted 

to A. Therefore, aigc^" is represented, modulo 0(/i), by -|- [k] + Since argc^" is itself 
a representative of [A] , the theorem is proved. □ 

If A = Aq is a whole Liouville torus, the principal terms in the quantization condition 
yield the well-known Maslov-Bohr-Sommerfeld quantization condition : "For all closed paths 7 
in Aq (or equivalently for the n standard generators of Hi{Aq) obtained through action-angle 
coordinates), 

iy"ao + y"^ + M7)|G27rZ + 0(/i)". (8) 
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5.5 Spectral parameter dependence 



The motivation of this paragraph is to malce the above results suitabfe for studying spectral 
problems, namely for investigating the solutions of the microlocal system : 

(Pl(/l) - ^i)^;, ~ 0, . . . , (Pnih) - En)uh ~ 0, 

according to the appropriate definition of microlocal equality (section 

In order to do this, the operators Pj in the quantization conditions have to be replaced by 
Pj — Ej. In the following analysis, we will need to state the problem in a slightly more general 
setting, where we replace Pj by an operator depending on E in essentially the same way as 

Definition 5.1 Let Pi{h) he a semi- classical integrable system on an n- dimensional manifold 
X, with principal symbols pi. A system Pf{h) depending smoothly on the parameter E G M" 
will be called regular deformation of Pi{h) if for each E, {Pf{h), . . . , P^{h)) is a semi-classical 
integrable system, with principal symbols pf , and there exists a C°° function f^{x) = f{E,x), 
/ : M" X M" — > M"- such that /° = id and /(•, 0) is a local diffeomorphism of (M", 0), such that: 

(pf , • • • ,Pn) = f^iPl, ■ ■ ■ ,Pn)- 

For such a regular deformation, the new Lagrangean foliation is globally the same as the 
original one Ac, with different labels : Ac = A^^^^^. 

Now we assume that there is a tubular neighborhood of A such that for E in a sufficiently 
small compact neighborhood of 0, the intersections A^ = A^ n O are non-singular submanifolds 
of A^. Then it is easy to verify that every result of this section extends to uniform solutions 
of the perturbed system Pf'{h) on i7. In particular we obtain a perturbed cocycle [c'^(/i)] that 
smoothly depends on E. It implies that the terms in the asymptotic expansion in powers of 
h of the cocycle [X^{h)], equal to the argument of [c^{h)], are locally smooth but possibly 
multivalued functions of E. At the principal levels, this gives a parameter dependent version of 
the Bohr-Sommerfeld conditions (^) : "for any family of loops 7^ drawn on Aq', 

i /" ao+ [ ac^ + Mt"")^ e27rZ + 0(/i)". (9) 

h J^E J^E 2 

Remark 5.2. Let P^{h) be such a smooth deformation. Then one can look at the family 
obtained by letting E = ht, where t is a bounded parameter (let's assume for clarity in this 
remark that n = 1; everything applies to the many-dimensional case as well). Using that in 
every local coordinates, the Weyl symbols of Pf is a smooth function of E, one can see that the 
operator 

P\h) = P^\h) 

is still a classical pseudo-differential operator, with uniform estimates for bounded t, but its 
principal and sub-principal symbols are of course not the same as those of P^{h). More precisely, 
if we denote by r^ the sub-principal symbol of P^ , then, because of Taylor's formula : 

f{E, x) = /(O, x)+E^iO,x) + --- + ^^^(0, x) + 0(E^+i), (10) 
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which holds uniformly for bounded E, the principal and sub-principal symbols and f* of P^{h) 
satisfy : 

p' = p' = f{p), and f* = rO + t^ (p). (11) 

OE \E=0 



Now, let us see what Bohr-Sommerfeld conditions (^) become. First of all, (|T^) shows that 
the principal Lagrangean Aq of the new system is independent of t and equal to Aq. Let us 
call «;* the sub-principal form on Ag. On the other hand, the family of Lagrangean manifolds 
A^ = {f^{p) = 0}, with E = ht, depends smoothly on h; let us call 5k* its deformation 1-form 
at /i = (see ^.2.2 ), defined by : 

dK^^.XfO/n) = —t—— (p). 
^ dE \E=Q ' 

Then paragraph 3.2.4 together with ([TT|) shows that, on Ajj, 

= + 5kK 

Now, if 7^ is a smooth family of loops drawn on A^, lemma implies that, when E = ht, 

/ ao= ao + h 5K^ + 0{h^). 

J'yE J.yO J^O 

Since it is clear that J^e = f^o + 0{h), we get : 

i / ao+ [ = \ f ao+ [ R* + 0{h). 

J-yE J^E n J^O J^o 

In other words, we exactly recover the Bohr-Sommerfeld quantization conditions we would have 
expected to find for the t-dependent system P*(/i). 

Note that this remark will mostly be used with f{E, x) = x — E, so that = — t. A 



6 Bohr-Sommerfeld for a focus-focus singularity 

6.1 Integrable systems with a non-degenerate singularity 

Let M be any symplectic manifold of dimension 2n, and /i, •••,/« a completely integrable 
system on M. The system is said to be singular at a point m G M if m is a critical point for 
the momentum map F. We shall always suppose in this work that m is of maximal corank as a 
critical point, which means that each function /j is critical at m. Without loss of generality, we 
shall also assume that for all i = 1, . . . , n, fi{m) = 0. 

The map Ti that assigns to a smooth function its Hessian at m is a Lie algebra homomorphism 
with respect to the Poisson bracket from the Lie algebra of smooth functions that are critical 
at m to the Lie algebra of quadratic forms on the tangent space TmM. Fixing any local set 
of Darboux coordinates, the latter is identified with the space Q(2n) of quadratic forms on 

M2- = {(x,0}. 

Every singular completely integrable system thus gives rise to a linear sub-algebra Ci;' of 
Q(2n), namely the sub-algebra generated by {^{{fi), . . . ,TC{fn)}- Note that Cp is always 
Abelian. 
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Definition 6.1 (Eliasson |23|) A singular completely integrable system with momentum map 



F admits m as a non-degenerate singularity if Cp is a Cartan sub-algebra ofQ{2n). 

Recall that a Cartan subalgebra of a semi-simple Lie algebra a is a commutative subalgebra 
c that is equal to its commutant and such that, for any element H £ c, adn is a semi-simple 
endomorphism of a. 



From the work of Williamson [43|, real Cartan sub-algebras of Q{2n) are classified according 
to the following scheme : 

Theorem 6.1 (Williamson) Let t be a real Cartan sub-algebra of Q{2n). Then there exists 
a symplectic linear change of coordinates in M^" and a basis gi, . . . , g„ of c such that each qi is 
one of the following : 

• qi = + ( elliptic type ) 

* qi = Xi^i (hyperbolic type) 

^ qi = Xi^i + Xi+iii+i 

Qi+i = Xiii+i — Xi+i^i 



(focus-focus type) 



If M has dimension 4, as we shall assume from now on, only four combinations are possible. 
We will restrict ourselves to the case of a focus-focus singularity. Note that this is the only 
case where m is isolated amongst the set of critical points of F. In particular, the singular 
foliation Ac = F~^{c) has, for small c, only one singular leaf Aq, and if J7 is a sufficiently small 
neighborhood of m in M, then O \ {m} is foliated by smooth Lagrangean submanifolds of M. If 
Aq is connected, this also implies that the critical value lies in the interior of the image of the 
momentum map F. The aim of the following paragraphs is to give a precise description of the 
singular fibration on a full neighborhood of Aq, generalizing the usual Arnold-Liouville theory 



for regular fibers. The results are not new, since they are quoted in |35] and proposition |6.2| can 
in principle be recovered from [^], but I believe that the following presentation gives a clearer 
account of what is needed for our purposes. 



6.2 Linear focus- focus 

We fix here the notations used when referring to a linear focus- focus system on = {(x, y, ^, r/)}. 
We shall always use Williamson coordinates, in which the system is generated by the quadratic 
forms qi = x^ + yr] and q2 = xrj — yS^. The Hamiltonian flows are denoted by and It will 
be convenient to use polar coordinates {r,6) for {x,y) and {p,a) for (^,7?), in which Xi and X2 
take a simpler form (see figure ||) . Of interest are also the complex representations zi = x -\- iy 
and Z2 = i + ir]. The flows of Xi and X2 are then respectively 

{t,{zi,Z2)) ^ {e^zi,e~^Z2), and 

{s,{zi,Z2))^{e''zi,e'''z2). (12) 

Similarly, we can identify the momentum map F with the complex- valued function F = qi -\-iq2, 
which gives 

F{zi,Z2) = Z1Z2. 

The singular leaf Aq = F~^{0) is the union of the coordinate planes Ps = {zi = 0} and 
Pu = {z2 = 0}, which are respectively the stable and unstable manifolds for the ^Yi-flow. 
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Figure 1: the linearized flow 



6.3 Geometry of the singular Lagrangean 



Singular points of focus-focus type are isolated, but it may occur that the singular Lagrangean 
Aq contains several of them. Since Aq is assumed to be compact, there will only be a finite 
number of them. However, we will always assume that Aq carries only one critical point, which 



is certainly a "generic" assumption (in |35], it is claimed that one can always reduce the number 
of focus-focus points by a small perturbation of the momentum map). This point will always be 
denoted by m. 

Recall that the usual theory of Arnold-Liouville gives action-angle coordinates in a neighbor- 
hood of any connected compact regular leaf. Here, we have compactness, but without smooth- 
ness. Another viewpoint would be to consider the non-compact punctured leaf Aq \ {m}; it is 
smooth, and invariant under the action of the joint Hamiltonian flow of the system (this holds 
because the critical point {m} is a fixed point). But there is no hope for the neighboring leaves 
to be diffeomorphic to Aq \ {m}, since they are compact... 

Nevertheless, the local structure of Aq near the singularity turns out to be sufficient data for 
its global description. 



Proposition 6.2 Let Aq be a singular leaf of the momentum map, carrying a unique critical 
point m of focus-focus type. Then the connected component of m in Aq is the image of a 
Lagrangean immersion of a 2-sphere with a double point. Deprived of m, it is an orbit of the 
Hamiltonian action of the system, with the structure of an affine infinite cylinder. 



Proof. First of all, note that we can suppose that Aq is connected. Indeed, this won't lead to 
any loss of information, since any other connected components are necessarily regular tori, for 
which the usual theory applies. 
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On Aq \ {m}, the action is locally free, hence by standard arguments each connected com- 
ponent of Aq \ {m} (there are at most two of them) is an orbit, on which all isotropy subgroups 
/ are conjugated. Therefore, this orbit is diffeomorphic to M?/I, which can be either M^, or the 
infinite cylinder M^/Z ~ 5^ x M (it cannot be a torus because it is not compact). 

The local structure of the flow makes up for this ambiguity. Indeed, according to Eliasson's 
theorem , the action is symplectically linearizable near m : there exists a symplectic chart in 
which the Hamiltonian vector fields of /i and /2 are linear combinations of the standard focus- 
focus vector fields; moreover, the involved coefficients form an invertible 2x2 matrix Mc which 
is locally constant along each fiber Ac. Now, the focus-focus vector fields are by assumption 
associated to the Hamiltonians qi and q2- Since q2 has periodic orbits in any neighborhood of m, 
the isotropy subgroup I is necessarily isomorphic to Z. The connected components of Aq \ {m} 
are therefore infinite cylinders. 

On such a cylinder are globally defined two infinitesimal generators of the action, Xi and X2, 
by the requirement that they are obtained from the initial generators Xj-^ and Xf^ by applying 
the matrix M~^. In Eliasson's coordinates, these vector fields are of course the standard focus- 
focus generators. The first one Xi is always transversal to X2, and hence describes an "axis" of 
the cylinder. Xi and X2 define the affine structure of the cylinder. 

Let us show now that Aq \ {m} is necessarily connected. From the connectedness of Aq and 
the local structure at the critical point m, we a priori know that Aq \ {m} has at most two 
connected components, which correspond respectively to the stable and unstable manifolds of 
Xi. Let X be a point on the unstable manifold, contained in a small S^-invariant neighborhood 
U of m. We let the fiow of Xi act on x so that when time increases, the image x(t) of x goes 
out of U. Now, we know from the dynamic on the infinite cylinder that x{t) must go out of any 
compact subset in finite time. Since the manifold with boundary Aq \ ?7 is compact, the image 
of X must necessarily enter U again. But it is clear from the local structure of the fiow that x{t) 
can approach m only via the stable manifold. Therefore, the stable and the unstable manifold 
are connected to each other, and hence they are equal. 

Thus, Aq is homeomorphic to a cylinder whose extremities are compactified in a unique point 
m, or a "pinched torus" , or a sphere with two points identified. 





























t — 







m 



Figure 2: topology of Aq 



Remark 6.1. We can readily deduce that 7ri(Ao) is isomorphic to Z, the integer coefficient 
counting the winding number along the "big axis of the pinched torus" . A 
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From the difFerentiable viewpoint, we see in the local form that Aq H ?7 is the union of 
two smooth submanifolds diffeomorphic to open discs and transversally intersecting at m. 
Therefore Ag is the immersion of a smooth 2-manifold with a unique double point. Since this 
manifold can be obtained by gluing a at each extremity of the cylinder Aq \ U , it is a sphere. 
□ 



6.4 Monodromy of the Lagrangean fibration around Aq 

The first description of the monodromy invariant for singular fibrations was carried out in 
llgj] , where the example of the spherical pendulum was studied. Later, Maorong Zou in |44| 



and Nguyen Tien Zung in |35| observed that the monodromy at stake was typical of a focus- 
focus singularity. In |35|, the computation of the monodromy is based on a claimed normal 
form for the foliation. We present here another proof of the result, relying on the local but 
Hamiltonian normal form of Eliasson. The constructions explained in the proof will be used for 



the quantization conditions (sections |6^ and 6.7). 



At the same time as I was writing this article, Duistermaat and Cushman were able to 
generalize the monodromy for non-Hamiltonian systems []l6| . 

Let us first briefly recall the definition, from a topological viewpoint. For any completely 
integrable system in dimension 2n, a regular value c of the momentum map F gives rise to 
a Liouville torus Ac, and thus to a lattice Lc = -ffi(Ac,Z). Because F is a smooth fibration 
near F~^{c), the union of all Lc is a locally trivial smooth bundle L over the open set Dj. of 
regular values of F. Fixing a point cq and a smooth basis of Hi{h.c) — Z", a local section of 
L near cq takes values in Z", and thus is locally constant. We therefore have the notion of 
parallel transport above a path in D^- The holonomy associated with it is a homomorphism 
from 7ri(L'r,co) to Aut(LcQ) — GL(n, Z), called the monodromy of L. 

Back to our focus-focus system, we take D^. to be a small ball around the origin in ~ C, 
minus the origin. Because its vri is Z, the problem is reduced to determining the monodromy of 
a simple loop around 0. 

Recall from last paragraph that for any focus- focus system (^1,^2)5 there is a local normal 
form (gi, 92) = f{pi,P2) near m, where / is a local diffeomorphism of (M^, 0), and g^'s are defined 



in symplectic coordinates around m by paragraph 6.2. This yields two particular vector fields 
Xi and defined on a full neighborhood of Aq by : 

{Xi,X2) = df.{Xp^,Xp^). 

The flow of X2 on A^ is, as we already saw, periodic of period 27r, so the orbit of any point x 
is a simple loop F^; on Ac depending smoothly on x. Its homology class 72(00) is therefore an 
invariant of the monodromy operator. Let us denote by 71 (cq) the homology class of a path 
starting on F^; in the direction of Xi, going back to F^, and closing up on F^; in the direction of 
X2. 

(7i(co)i 72(co)) is a basis of Lcq depending locally smoothly on cq; the monodromy of the 
system is expressed on this basis by the following proposition : 



Proposition 6.3 ( [|35f| ) Let (^1,^2) a completely integrable system with a singular leaf car- 
rying a unique singular point of focus-focus type. Then for any cq close enough the the critical 
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value G M^, the monodromy matrix, expressed in the basis (7i(co), 72(00)) of Lcq, is equal to 




Here e is the sign ofdetM, where M G GL(2,M) is the unique matrix such that : 

{n{qi),n{q2)) = M.{n{pi),n{p2))- 

Proof. The idea of the proof, aheady apparent in the construction of the basis (71(00), 72(00)), 
is to show that all contributions to the monodromy can be concentrated in a small neighborhood 
of the critical point m. 

Let U be such a small ball around m, taken to be 5^-invariant. The intersection of dU with 
Ao is composed of two circles and Tg respectively drawn on the local unstable and stable 
manifolds at m. Let us take a point x on F^, and let S be a small section transversal to Ao at x. 
The tangent space to Ao at x is exactly the kernel of dF. Thus, the momentum map F realizes 
a local diffeomorphism from (S,x) onto (M^,0). We can then assume that S is a preimage of a 
small open disc D C Dr around 0. 

Let Q be the invariant open subset of M obtained as the union of orbits of points in E. 
It is a neighborhood of Aq, composed of the union of Ao and invariant tori. Because of the 
compactness oi 0,\U, any orbit of some point y E S close to x intersect dU in two circles close 
to F„ and F^. The restriction of F to i7 \ [/ is therefore a trivial fibration over D whose fibers 
are finite length 5^-invariant affine cylinders, where the affine structure that defines the "axis" 
of each cylinder is given by the vector field : 

n\U ^Dx ([0, 1] X S^) 

F (13) 
D 

For instance, the circles F^ and F^ are represented in this trivialization by (0, (0,S^)) and 
(0, (1, S^)) respectively (see figure ^. For any y £ Q\U corresponding to the point (c, {t, r)) in 
the above trivialization, we denote by Us] the A'l-path corresponding to the axis (c, ([0, 1], t)). 
As we can assume that U lies inside an open neighborhood of m where Eliasson's normal form 
applies, yu and yg have complex coordinates {yu/s,iiyu/s,2) satisfying 

¥^^yu/s,2 = /(c)- 

Note that the matrix M in the statement of the proposition is equal to df{m). 

Remark 6.2. We don't actually need the full statement of Eliasson's normal form, because 
higher order terms in the normal form wouldn't be of any trouble for the local computation of 
the monodromy (see [^). A 

If c / 0, using formula (|l2|) for the fiow, ys can be joined to yu by letting the flow of qi act 
during a time t{y) = ln{\yu^i\/\ys,i\) and that of q2 during a time 

s{y) = aig{yu,i/ys,i)- 
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Figure 3: The fibration F on Q\U 



By this construction, we get a locally smooth family of loops on Ac, whose homology class at 
Co, 7i(co), together with 72(00), form the basis referred to in the statement of the proposition. 

Let be the axis through x (by definition, x„ = x is on the local unstable manifold 

Pu and Xs on the local stable manifold Ps). If y tends to x, then y„ and Us tend to Xu and 
Xs respectively. Therefore, since Xu,2 = Xg^i = 0, s{y) ~ aT:g{xu,ixJ^) + arg(/(c)) as y tends 
towards x. 

Now let C be a small loop around x in S such that its image by F is a circle in D through cq, 
oriented in the trigonometric sense : letting y run around C, and with it, the loop 72(c), amounts 
to increasing arg(c). Therefore, after such a revolution, s{y) get increased or decreased by 2tt, 
depending on the diffeomorphism / being respectively orientation preserving or reversing. This 
means that the loop obtained is equal to the initial one composed by respectively plus or minus 
an orbit, which proves the proposition. □ 



6.5 Microlocal focus-focus 

We are now given a self-adjoint semi-classical completely integrable system Pi{h), P2{h) whose 
associated classical system presents a focus-focus singularity at a point m, as defined in the 
preceding section. 

The aim of this section is the complete analysis of microlocal solutions around the critical 
point m of the system : 

in which P^ is a regular deformation of Pj with a parameter E varying in a compact neighbor- 
hood of e (see definition [ST]) . 



The problem of globalizing these solution to a neighborhood of the whole critical leaf Aq will 
be dealt with in the next section. 

For now, the crucial point is the following : 
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Proposition 6.4 The space of microlocal solution of (\l4 ) at m has dimension ^1. In other 
words, if , {h,E) £ T is a non-trivial solution of ^1^ ) at m, then for any other such solution 
, there exists a constant C^{h) G C^, depending smoothly on E, such that : 

vj^ ~ C{hfuf at m. 

To prove this, we use a parameter dependent version of the normal form of |^9[ : let Qj{h) be 
the standard focus-focus operators on R^, that is, written in polar coordinates {r,9) : 

^ h , d 1 ^ h d 

(5i = -(rTr + l) and Q2 = - — . 
I or I Oil 

For notation convenience, we will write P{h) for the vector {Pi(h), P2{h)), and similarly for 
other quantities. Let p be the momentum map for the system (pi,P2)- The target space M? of p 
will be identified with C by p = + ip2- As we already saw, in suitable symplectic coordinates 
near m, there is a local diffeomorphism F of such that p = F{q). Thus, 

p^ = f^oF{q). 

The semi-classical normal form that we get is the following : 

Lemma 6.5 ([^9|]) There exists a unitary Fourier integral operator U^{h), a microlocally in- 
vertible 2x2 matrix M.^{h) of pseudo- differential operators, and e^{h) S Ch, everything de- 
pending smoothly on E, such that, microlocally on a neighborhood of m : 

admits a semi-classical asymptotic expansion of the form : 

e^{h)^Y.^^el 

fc^o 

The first two terms satisfy : 

e^ = F-'o{f^)-\0), 



e? = -M"\(r(m)). 

• r(m) is the value at m of the sub-principal symbol of P. 

• M is the value at of principal symbol of Ai^. It is also defined by M = dF{0), or, equiva- 
lently : 

d^p^{m) = M.d^q{0). 



Remark 6.3. We will show later (paragraph |6.7.4 ) that the expansion e|'(/i) = 3f(e^(/i)) is, up 



to 0{h°°), independent on the choice of a Fourier integral operator U{h). For the other "semi- 
classical invariant" ^^{h) = 9^(e^(/i)), this might not hold; however, will will prove thanks to 
corollary |6.1C that its formal infinite germ at = is indeed invariant by a change of Fourier 



integral operator U{h). A 
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Remark 6.4. If we apply this lemma to E = ht for a bounded n-dimensional parameter t, we 
can apply remark |5.2| and express everything in terms of t. We obtain an asymptotic expansion 

e*(/i) =^ e^^{h) satisfying : 

e* = 0, ~e{ = -M'\ (r^{m) + ^ (0)i 



dE \E=o 

A 

Thanks to this lemma, the solutions of (|l^ get transformed into those of simple differential 
operators on for which we have explicit solutions. We can now turn to the proof of proposition 
6.4. Assume that we know an exact solution of the system uf^^^^ G 5'(M^) (this will always occur 
as soon as the system admits a microlocal solution - see proposition ^^). Let uj^ be any other 
microlocal solution near 0. 

Let Tr be the "vertical strip" {r < R} in the cotangent bundle r*M^ = {{r,6), (p, a)). The 
Lagrangean leaves of the system in T^, deprived of a small ball B around 0, are non-singular, 



hence from proposition there exists a constant C (h) such that microlocally on T/j \ B, 

(remember that such a formula involves uniformity with respect to E\ 

As in [0], we shall use this exact solution to construct a local (mod 0{h°°)) solution, 
microlocally equal to uj^ near the origin. Let us define by : 

^hv^ = xThuf + (1 - x)-^hC^^xf_„ 

where x is a compactly supported C°° function with value 1 in a neighborhood of 0. Then 
and are microlocally equal near the origin. Furthermore, is a microlocal solution of the 
system at each point of Tr. 

Now there exists differential operators Pfiji) such that for any w G 5'(M^), ThiP^vS) = 
P^J^hW. Differential operators preserve the support of distributions; therefore, since ThV^ = 
•^^^(fact f^'^ from the origin, we have there P^J^hV^ = ^hiPfuf^^^J = 0. So we have 

Pfv^ ~ 

near each point of Tr, and 

:F,{pfv^) = 



far from 0. We deduce from this that Pfvj^ satisfies the hypothesis of lemma 4^ on T^j, that 
is, 

Pfv^ = 0{h°^) on TKn{p = 0}. 

We have thus at our disposal a solution of the system Pfvj^ = w^^ on with w^^ = 
0{h°°) for r < R. On \ {0} this is a regular linear differential system of the first order, whose 
solutions are of the form 

= Cste'^{h).ut^, + (15) 
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^^act is a solution of the homogeneous system, and w^rticuiar is a particular solution of the system 
that can be computed from u^^^^ by the method of "variation of the constant" . Explicitly, if we 
write 



we find that the function must satisfy the following system, expressed in polar coordinates : 



dg^ _ 


i 




dr 


h 


''exact 


dg^ _ 


i 




~5e- - 


h 


"exact 



If (ro, ^o) is any point with < ro < R, is integrated by the formula : 

g^{r, 6) = C a^iro + t{r - tq), + t{d - 0o)).(r, d)dt + ^^(ro, ^o)- 
Jo 



It is easy to check from the explicit formula of uf^^^ (see proposition |6.6| ) that is uniformly 
0{h°°) in the disc r < i?; so with Cste^{h) = g^{r(), Oq) we get in (15) that w^^^j^^i^^ is uniformly 
0{h°^) near the origin. 

The distribution zj^ = vf — C ste^ .uf^^^^ ~ ''^^rticuiar is a solution of Pfz^ = w^^ with support 
concentrated at the origin. Therefore, it is a finite sum of derivatives of the Dirac distribution. 
Now, the identity Pi'z^ = 0{h°°) readily shows that the the coefficients involved in zj^ are 
uniformly 0{h°°), and so z^ is. 

Finally we are left with vj^ = Cste^.uf^^^ + 0{h°°) near 0, which implies vj^ ~ Cste^.uf^^^^ 



near 0. This, provided we prove proposition |6.6|, achieves the proof of proposition 3.4. □ 



Remark 6.5. The proof shows that the constants C^{h) and Cste^{h) are microlocally equal. 
Besides from the preceding arguments, this essentially comes from the local connectedness of 
the Lagrangean leaves Af for S / 0. A 

It remains now to check : 

Proposition 6.6 We consider the system (^Jj on with P^{h) = Qj — e^{h), as previously 
defined. The following conditions are equivalent : 

1. there exists a non-trivial microlocal solution , {h,E) G F o/ ( jl^ on a neighborhood of 
the origin; 

2. up to 0{h°°) modifications of the , there exists an exact solution u^^^j of ( [7^ in M^\{0} 
for {h,E)eT; 

3. £Z + 0{h°°) uniformly for {h, E) E F. 

Proof. Equivalence between 2. and 3. comes from the fact that at each point of \ {0} the 
solutions of the system are spanned by the function : 



-ie^^i^V^^ (16) 
r 
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Therefore, there exists a solution around the origin if and only if G Z. 

Since uf^^^^ as defined above is an admissible functional in depending smoothly on E, and 
a solution of the system on M^, it is also a microlocal solution. Hence 2. implies 1. 

Now let uj^ be a non-trivial microlocal solution on a neighborhood Q of 0. From the unique- 
ness of microlocal solutions at points different from the origin (proposition |5.1| ) , we deduce that 
on every simply connected open subset fio, of $7 \ {0}, there is a constant (C^)°(/i) such that 

Then, going around 0, the (C^)"'s must microlocally coincide. Hence the e*~^ must all mi- 
crolocally coincide, uniformly for (h, E) G F, which implies 3. □ 

Put together, the two propositions give the following result ( "first quantization condition" ) : 



Theorem 6.7 The system ( fl^j has a non-trivial microlocal solution (with {h,E) G on a 
neighborhood of m if and only if 

G Z + 0(/i°°), {h,E)eT. 

lb 

Any other solution is a multiple of it. 



Remark 6.6 

on the base space. Here, we can rewrite the half-density m, 



As we already mentioned, microlocal solutions should be viewed as half-densities 

Lct\dxAdy\^ as 



ln|rl 



dx A dy 



and for r ^ 0, its modulus 



dxAdy 



is exactly the pull down of the canonical invariant Liouville 



half-density p{X) on A^. In this sense. 



is a normalized classical oscillatory half-density on 



every open subset A^ of Aq" that projects diffeomorphically onto X (thus excluding, for E = 0, 
the vertical (^, r/)-space in T*X). Concerning the phase function (fE = ef o In |r| -|- efg^j a handy 
tool for its expression is the use of complex coordinates zi = x + iy and 22 = ^ + irj; we also 
write ej^ = ef q + ie^n, which leads to 



2,0' 



9^(eo Inzi) 



(Note that EqIuzi is not well-defined, but thanks to the assumption on e2'(/i), its real part is 
indeed well-defined modulo 27r/iZ.) We readily deduce that the graph of ip^ in T*X is the set 
of (zi, Z2) such that Z2 = 2^ = e^/zi. Thus Af = {(zi, 22) G T*X, Z1Z2 = ef }. Note that 
this also holds if -E = 0, for then ef = and Ag consists of the union of both coordinate planes 
zi = and Z2 = 0- A 
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6.6 Regularization of [k] 



Before stating our principal result in the next section (theorem |6.9| ), we need to describe one 
more geometrical ingredient, related to the divergence of the integral of the sub-principal form 
K along paths through m. 

We still let {Pi{h), P2{h)) be a semi-classical completely integrable system with a singular 
principal Lagrangean carrying a unique singular point m of focus-focus type. Here again, we 
identify the target space of the momentum map p with the complex plane C. Thus p = pi + ip2- 
For small c S C, we denote by the Lagrangean leaf = {p~^(c)}, and usual the 

sub-principal form on Ac. 

Let 7c be a smooth family of loops on Ac such that 70 is a simple loop through m starting 
on the local unstable manifold of (Ao,m), and coming back to m via the local stable manifold. 
We know that the non-triviality of the monodromy implies that the map c — > 7c cannot be 
univalued, but proposition gives a precise control over its multivaluedness. The purpose of 
this paragraph is to control the integral of the sub-principal form along such a family of paths, 
and in particular to determine its divergence principal part as c ^ 0. 

Let M be the 2x2 real matrix such that TCi^p) = M.7i{q), where q = qi + iq2 (for the 



linearized flow, we always use the notations of |6.2| ). In accordance with lemma 3.5 we let 
where r = ri + ir2 is the joint sub-principal symbol of the system. 

Proposition 6.8 For any multivalued smooth family 7c of loops on Ac such that 70 is a simple 
loop through m starting on the local unstable manifold of {AQ,m), and coming back to m via the 
local stable manifold, then 



the limit : 



lim (j Kc + lH(eiln(M-ic))j 

along a ray '^^ ^ 



exists, 



and its class modulo 27r3(ei)Z actually does not depend on 7c, provided 70 satisfies the 
hypothesis; 

finally, this class is also given by the formula : 

170(1^) = , / i^ + ^{ei)ln{rAoPBo) + 3{ei){OAo-aBo)] , (17) 

(S,t)^(0,0) yjAo=-yo(s) J 

where (r, 6, p, a) are the polar symplectic coordinates near the origin defined in paragraph 



6.i. 



I-yo^K) will be called the principal value of J^^ k. 
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Proof. Let ([71(c)], [72(c)]) be the smooth basis of 7ri(Ac) defined in proposition |6.3| . The 
homotopy class of 7c decomposes into 

[7c] = ni(c)[7i(c)] +722 (c) [72(c)], 

for integers ni(c) and n2(c). On Aq, 72(0) degenerates into a trivial loop, and by hypothesis we 
have 

[70] = [71(0)]. 



But we know from proposition 3.S that the group bundle 7ri(Ac)/(72(c)) ^ c is trivial near 
c = 0, with fibers isomorphic to Z. Therefore, ni(c) must be identically equal to 1 for small c. 

If c is constrained into a sector A of the plane, then the family 7c must be continuous at 
zero, which implies that n2(c) is constant for small c. Since the integral of Hc along the Af2-circle 
72(c) is equal to — 27rr2(m) + 0(c), we have 



lim / Kc = 27rg(ei). 

ci*A -^"2 (c) [72(c)] 

This, provided we prove the first point of the proposition, shows the second point. It remains 



to show the convergence of /-^^(c) '^c, and the formula 17. We shall still assume that c stays in a 
sector A and by 7c we denote a smooth representative of the first basis element [71(c)]. 

Let / be the function on (IR^,0) such that, in some neighborhood of m and some local 
symplectic coordinates around m, p = f{q), so that M = df{0). If we let c = /(c), the problem 
amounts to investigating the limit of 



/ Kc- + $n(eiln(c)), 



because ln(M ^ f{c)) = ln(c) + 0(c). Now, Kg is defined by 
which also reads, in Q : 

This means that Kc is, in fi, the sub-principal form on the c-level set of a system with principal 
symbol q and sub-principal symbol {df{q))~^.r (which we call r again). So let us work with 
these data and forget the tilde. 

As we shall see, we can even reduce the problem to a constant sub-principal symbol r(0)(= 
— ei). This is due to the following "critical Poincare lemma" (see [^]) : there exists a function 
(7 on a neighborhood of G T*1R^ and smooth functions hj on (]R^,0) such that 



Tj = dg.Xq^ + hj{q). 
This implies that on Q, for non-zero c, the function 

Kc = -g- hi{c) ln(r) - h2{c)9 ( mod 27r/i2(c)) 
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is a primitive of Kc on Ac = {rpe'^^"'~^^ = c}. Suppose now that Aq and Bq are points on 70 lying 
respectively on the local unstable and stable manifolds of Aq. Then there exist smooth families 
of points Ac and on QCi 7c([0, 1]) such that limc^o(^c) B^) = {Aq, Bq). We can write : 



[ Kc= [ Kc + KciAc) - KciBc 



Be 

Kc + g{Bc) - g{A^ - hx{c)ln(r aJtb^) - h2ic){0A, 



Ac 

which can also be written 

"B, 



/ Kc + g{Bc) - g{Ac) - hi{c)ln{rAcPBc) - h2{c){0Ac - as^ 
J Ac 



+ {hi{c)ln{\c\) -/12(c) arg(c)). 

Everything but the last parenthesis has a smooth behavior as c — > in A, and that last paren- 
thesis is nothing else than 9^((/ii(c) + «/i2(c)) In c). So there is a smooth function I{Ac, Be) such 
that : 



IiAc,Bc)= / Kc-9\{{hi{c) +ih2{c))lnc). 



Remark 6.7. Because of the monodromy matrix ^ | ^ ) (proposition |6.3| ), as c loops in 
the positive sense around the origin, J Kc gets increased by 



/ Kc = -27r /12(c) 
-'72(c) 



whereas 9^((/ii(c) + i/i2(c)) In c) = hi{c)ln{\c\) — /12(c) arg(c) gets increased by the same quantity. 
This shows that I{Ac, Be) is actually a univalued smooth function around c = 0. A 

Because (/ii(c) + i/i2(c)) = — ei + 0(c), we deduce that if we fix a determination of Inc and 
make c tend to in A, then 

9^((/ii(c) +i/i2(c))lnc+ eilnc) 

Therefore, I-y^iK) exists and is equal to I{Aq, Bq). Note that it is independent of Aq and Bq, and 
hence equal to the limit of I{Aq, Bq) as Aq and Bq tend to the critical point m. That means that 
the function g, being smooth at m, plays no role, justifying our assertion that the sub-principal 
symbols can be taken to be constant near m. This remark also yields the useful expression of 
I-yo{K) involving only objects living on Aq given by formula (0). □ 
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Figure 4: regularization of k 



6.7 Global solutions 

Here is finally the main result of our work, namely the condition under which there exists a 
microlocal solution of the system (|l^ on a neighborhood of the whole critical Lagrangean Aq. 
The hypothesis are the same as in paragraph |6.5| . From now on, n is a generic integer, and not 
the dimension of the ambient symplectic manifold, which is 4. 



Theorem 6.9 Let = ^f{h) be the constants defined by lemma |67 

• There exists a T-valued constant (A^)°"' = (A^)°"*(/i) admitting an asymptotic expansion in 
^ — 1 powers of h whose terms are smooth functions of E, such that the system ^1^) admits a 
microlocal solution u^, {h,E) G P, on a full neighborhood of Aq if and only if the following two 
conditions are satisfied : 

1. 'first quantization condition" : 

ilR^Z + Oih^); 

2. "second quantization condition" : 

(XT + "I - f H2h) - 2 argr (^Jl/h+l+^^ ^ 2vrZ + 0{h^), 

where -^-^ ~ n G Z. 

• Moreover, the first two terms in the asymptotic expansion of (A^)°"'(/i) are the following : 

• for non-zero E, 

+ ^^K^ + ef;i/n|c| - enlarge + /x(7f)|j , 
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where "yf is a smooth family of loops on Aq such that satisfies the hypothesis of propo- 
sition \6.^ , c = efg + ie2Q, and /x is the Maslov cocycle on Aq . For small E, /i(7f') is a 
constant integer. 

• forE = 0, 

where k = kP is the sub-principal form of the system (Pf , ^-'2 )• 

The last point of this theorem follows from the expression for non-zero E (recall that A'^ is then 
smooth) via proposition |6.8| : indeed, if we let c = (/■^)~^(0), we have = Aq', and the family 
of loops 7a satisfies the hypothesis of proposition |6.8|. The 1-form is defined by : 



{^^.X^o,^^.X^o) = {df^)-\{rf,ri). 

Writing ii^ as a function of c, this shows that is the sub-principal form of a system with 
principal symbol (^1,^2) sub-principal symbol ((i/^)~^.(rf', r|'). The value at m of this 
sub-principal symbol is r*^(m), hence proposition implies that 



/ K-^ + in(e?lnM-ic) 



has a limit I as E ^ 0. Because ef tends to and c = is tangent to M at c = 0, this 
limit is equal to that of 



/ K^ + fn(eflnc) 



Moreover, we saw in ( |l7| ) that it only depends on the value of the sub-principal form on Aq. 
But, on Aq, {df^)~^ .{ri,r2) = r^. Therefore, / is equal to the principal value of a system 

with principal symbol and sub-principal symbol r^. □ 



The rest of the proof of this theorem will be split into two propositions |6.11| and |6.14 . 

We know from the geometry of the Lagrangean fibration (section |6.4D that a neighborhood of 
the principal Lagrangean Aq also contains regular tori Aq', E (E small) for which we already 



derived the quantization conditions (theorem 5^). Therefore, one of our further tasks will be 



to show how these regular quantization conditions are contained in the "singular" quantization 
conditions announced above. This link is provided in paragraph |6.7.4 



Finally, this theorem can be applied to "zoom" to a 0{h) scale : let, as in remark 5.2 , 
{Pi{h), P2{h)) be the semi-classical completely integrable system defined by 

Pjih)':^'p^\h), 

where i is a 2-dimensional bounded parameter. From the constant e^{h) = (ef'(/i), e|'(/i)) given 



by lemma 3^ we define e*(/i) *== j-e^^{h). This e* admits an asymptotic expansion in positive 



powers of h whose coefficients are smooth functions of t. The principal term is given by : 

~t _ ^ , , 
dE \E=o 
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that is to say, in view of (pA|), 

el = -M-^.f\m). 

Here f* denotes the joint sub-principal symbol of the system (PfjPg)- Recall that its principal 
symbol \s p = and does not depend on t. 

Similarly, the holonomy coefficient (A^)°"*(/i) ~ Y^^=-i{^^)'k^ ^ t-dependent version 
(A*)°"'(/i) whose first terms are : 

We know from the theorem that ^(A'^)'^^{ = J^o oq- For the other term, using the identification 
= ef + , have, for non-zero E : 

a(A^)°_"{ d f E d 



dE 



Therefore, 



O\out 



with E = ht. Using lemma 3.2, one easily shows that, for fixed /i = /iq 7^ 0, J e is 



equal to f b (5k^, where the 1-form 6k^ is given by 

f)fE 



E 



Therefore, the sum 5k + k is the sub-principal form of a system with principal symbol p 
and sub-principal symbol ^^-t + r^, whose value at m is equal to f^{m). It follows that 

hm + K^)+m ((ef + ^.t) In 60^) = I^oiSn' + k). 

But + k is the 1-form on Aq that takes the value ^^=0'^ ~^ '''^ °^ ('^p!j' ' '^pO ) : it is exactly 
the sub-principal form of the system (Pf, Pg)- We have proved the following result : 
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Corollary 6.10 The system (P^jPg) defined above admits a microlocal solution on a neigh- 
borhood of the critical Lagrangean Aq if and only if the following two conditions are satisfied : 

1. "first quantization condition" : 

4(/i) G Z + 0(/i°°); 

2. "second quantization condition" : 

(A*)-* + n| - e\ H2h) - 2 argT (^A±1±J1^ ^ 27rZ + 0(/i~), 
where e^{h) n £ 7^. 

• Moreover, the first two terms in the asymptotic expansion of (A*)°"'(/i) are the following : 

^^^ao+ /^o(K*) + /i(7°)^, 

where is a loop on Aq that satisfies the hypothesis of proposition \6.^ and R is the sub-principal 
form of the system. 

Remark 6.8. This corollary makes it easy to show that, provided the first quantization condi- 
tion is satisfied, then the "semi-classical invariant" e\{h) is intrinsically defined by the microlocal 
behavior of the system around the critical point. This gives an a posteriori justification of the 



claim following lemma S.5 



Indeed, let n be the integer such that eK/i) ~ n (we shall prove in paragraph |6.7.4 that e2(/i) 
and even ef'(/i) - are intrinsically defined by the system around m, hence so is n). 

For bounded t, e{{h) is uniformly bounded as h ^ 0. As will be shown in the proof of 



proposition |6.11 , for any integer n G Z, the function p^i^^^"^]t|?^|^ analytic for x £ 



> 2 \ 

which implies that the function 2argr(|(ix + 1 + n)) is locally analytic. Therefore, plugging 
in X = e\{h) we obtain a semi-classical expansion in non-negative powers of h, whose terms are 
smooth T-valued functions of t. 

This in turn implies that the value of the cocycle [A*(/i)] on which gives the "second 
quantization condition" of corollary |6.10| , can be written under the following form : 

[A*(M](7?) = C\h) - e\{h)lnh + 0{h^), 

where C^{h) is a semi-classical expansion in ^ —1 powers of h whose terms are smooth functions 
oft. Therefore, —e\{h) is uniquely defined as the coefficient of In /i in the expansion of [A* (/i)] (7]^). 

Now, if the system is changed outside a neighborhood of the critical point, then only (A*)°"' 
is perturbed, which means only a semi-classical perturbation that has an asymptotic expansion 
in powers of h, which does not affect the Ink coefficient. A 

Remark 6.9. In dimension 1, the singular Bohr-Sommerfeld conditions of Colin de Verdiere 



and Parisse has been recently extended [13| to handle the case of multiple critical points of saddle 
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type on the critical Lagrangean. In our case as well, if we restrict to focus-focus singularities, I 
believe that the results presented here can be extended to treat k > 1 focus-focus points on Aq. 
This basically holds because the topology of such a Aq is just a /c-times pinched torus, and the 
monodromy if equal to the k-th power of the simple monodromy. However, including other kinds 
of singularities should require a more difficult analysis, starting from the topological description 
of Aq. Understanding the work of Fomenko (e.g. [25|), should probably help coping with this 
difficulty. A 



6.7.1 The microlocal bundle 

As a starting point, we shall assume that the system fulfills the first quantization condition 
(theorem 6.7). This is of course a necessary condition for the existence of a global solution. 

From this assumption, there exists a microlocal solution (u^)^ on a neighborhood of the 



critical point m. As in paragraph ^.2| , one can thus form a microlocal line bundle £h{^o ) over 
the topological manifold A^, with well-defined holonomy, by means of the trivializations (0). 
Of course this yields a family of line bundles depending on the 2-dimensional parameter E, but 
gathering the results of paragraph |5.5| and proposition |6.4| we see that the trivialization functions 
in can be chosen to depend smoothly on E. 

As a consequence, the holonomy of £h{^o) depends smoothly on E, in the sense that the 
elements (c^)"^ of the Cech cocycle c^{h) representing it are C°° functions of E with values in 

From this picture, we already know the abstract form of the quantization condition we are 



looking for : it is exactly the parameter dependent version of proposition 5.5. In other words, 
there exists a microlocal solution uj^ on a neighborhood of Aq if and only if for any family 7^ 
of loops in A^ that lie in that neighborhood, the value of c^{h) on 7^ is uniformly 1 + 0{h°°) 
(here, we are viewing as a homomorphism from 7ri(A^) to C^*). Let us now explicit this 
condition. 



Let i7 be a neighborhood of m in T*X in which the normal form of lemma |6^ applies. Then, 
as the system admits a microlocal solution on Q., any loop contained in iln Aq' has microlocally 
trivial holonomy. But for every small E ^ {), Q. contains a generator ^2 of '''"1(^0')' namely the 
one associated with the rotational symmetry of the fibration (it is an integral curve of the 
vector field X2 of paragraph [6.31 ). Such a ^2 ^-^so exists for S = but it is then homotopic to 
{m}. Anyhow, this yields : 

c^(7f) = l + 0(/i°°), 

uniformly for small E (this is of course due to the fulfillment of the first quantization condition) . 

Thus c^ restricts to the quotient group 7ri(A^)/(7^). For fixed every 7ri(A^)/(7^) is 
isomorphic to Z. What's more, proposition |6.3| implies that the group bundle 7ri(AQ')/(7|') — > 
U B E over a neighborhood U of in M2 is trivial : 

Let T'f' be the generator of 7ri(A^)/(7^) equal to {E, 1) in the above trivialization. The "second 
quantization condition" is now the requirement that c^(7'f') = 1 + 0{h°°). 
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Let us now examine that statement more closely. Recall that the line bundle Ch{A.Q) was 
defined through a set of local trivializations subordinated to a cover UQq of a neighborhood of 
Aq in T*X, such that on every existed a non-trivial microlocal solution (u^)" of the system. 
We can assume that the critical point m is contained in a unique element of that cover, which 
we will denote by Qq. Let us pick up a simple loop 7i(t), t £ [0,1], in Aq representing the 
quotient class and enumerate Qq, Qi, . . . ,r2£,Qo the O^'s encountered by 7i(t) on its way 
from m back to m again. We will here choose so that it goes away from m via the local stable 
manifold of Aq and returns back to m via the local unstable manifold. Also, we will suppose 
that U ill U C ^■ 

Now for every fixed E, we choose a representative 7f' of 7f'. Because of the triviality of the 
fibration outside m (13), the paths can be homotopically deformed so as to be close to j^, 



in the sense that their images enter IJj=o 

Then by definition, the value of on is the product : 

c^(^f ) = (c^)0'i(c^)i'2 . . . ^^Ey-i/^^Eyfi^ 

We are going to calculate this expression by splitting it in two parts : a local one - the product 
(c^)°'^(c^)^''^ - and an outer one - the rest. For this purpose, we are free to fix the solutions 
{uj^y on used for the local trivialization of £/j(Aq'); let us make the following choice : 

• choice 1 : for i 0, 1,£, {uj^Y is any classical oscillatory integral constructed as in proposition 
^.7| , with invariant half density given by the canonical Liouville density on Aq'. 

For i = 0, 1,£, we make use of the local analysis of the preceding paragraph To begin 



with, we transform the system on J7 to the normal form of lemma 3.5 by a unitary Fourier 
integral operator U{h). The microlocal solutions of the transformed system on 0, are spanned 
by: 



1 

E ?-J-lnr 7-2-6 

exact 



r 

Because of the first quantization condition, we know to be microlocally equal to an integer 
n S Z. For any real e, let u^^n be the tempered distribution on : 

_ \ ine ielrvr 

r 

so that u^^^^^ = u^E ,E . The distribution n is a solution of the system 



"1 '-2 
h ' h 



Id _ 

j Q0'^£,n — ^^e,n 

i(r|: + lK,. 



(18) 



h '' h h '' h 

system on $7. This allows us to complete our choice of local solutions by the following : 



on M^, and it is easy to see that its Fourier transform is a solution of the same system with 
£ changed to —e. Therefore u^e ^e and !F'i^'^{u_^e ^e) are two solutions of our transformed 

s all 

[ «)0~C/(n^^) onl^o 

(uf)i ~C/(n^^) on 1^1 

(nf)^~[/(.FfC_^^) onQ, 



choice 2 : < 



h ' h 
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6.7.2 The local holonomy at m 



We give here an explicit expression for the product (c^)^''^(c^)^'^, which is defined to be the 
constant C^{h) such that, on Qq n O^, («f )° ~ C^{u^y. 

Proposition 6.11 The unique microlocal constant C^{h) such that (uj^)^ ~ C^{uj^Y on ^IqCi^Ii 
is given by : 

^ I'ief/h + l + n 



-ief /h + 1 + n 



with ~ n S Z. 

Proof. Applying the Fourier integral operator Uh, we see that is also defined to be the 
constant such that, microlocally on VLq n , 



h '' h h h 



It turns out that the constant thus defined is a microlocal version of an exact problem 
concerning the homogeneous distributions on M^. 

We know that the tempered distributions u^^n are solutions of the system ( |T8|) in R^. They 
are in fact the only ones - up to multiplicative constants, naturally. To see this, one first restricts 
to \ {0}, where the result is standard, and then prove that the system (18) cannot admit 
linear combinations of derivatives of the Dirac distribution as solutions. Now, „ is also 

a tempered distribution and a solution of the same system, which leads to the first point of the 
following proposition : 



Proposition 6.12 ([Q, ^g]) 1. For any (e,n) G MxZ, there exists a unique constant C{e,n) 
such that : 

2. C{e,n) has the following expression : 

C(e, n) = rl"l2'^ ^ . 2,^, ■ (19) 

^ ' ' -ie+l+\n\ -J ^ ' 



The computation of C(e, n) is done by testing u^^n on derivatives of Gaussian functions (^)"e~ 

cfz) 



and (J=)"e^2' {^z = re^d). 



Remark 6.10. Apparently, this formula appeared for the first time in Tate's thesis, during 
the year 1950. However, it has not been published until 1967, when it appeared in the book ||^, 
chapter XV]. But at the same time as Tate wrote his thesis, the idea was in the air. Gelfand was 
studying homogeneous distributions in the real and published a closely related formula in [27|. 
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The complex version, that is, exactly the result mentioned in the above proposition, appeared 



in the addendum of the French edition of the book by Gelfand & al. on distributions [26|, 
where it is claimed to be published for the first time; the author was certainly unaware of Tate's 
work. Also interestingly enough, few years before (1951), a slightly different version of the same 
formula can be found in Bochner's work |2|. It is actually a generalization to homogeneous 
distributions on M'^. A 

Now, let us (momentarily) denote by D{e,n) the quantity : 



D{e,n) = r"2^" ^ ■ 

^ ' ' —i£+l+n 

The Gamma functions involved are analytic functions of e that, for negative odd n, have simple 
poles at e = 0. The quotient D{e,n) is therefore analytic - with no pole - on the real line (we 
won't consider here analytic continuation to complex e). More precisely, for n = — 2m — 1, m ^ 0, 
as the residues of the numerator and the denominator are the same, we have : D{0, —2m — 1) = 
— For non-zero e, we can use the complement relation to prove that D{e,n) = D{e, —n). 
Then of course this also holds for e = 0. Since D{e, n) and C(e, n) coincide for n ^ 0, we deduce 
that they are equal for all n G Z. 

Finally, using the relation J^hu{(,,r]) = jiJ^u{j^, ji^), we terminate the proof of proposition 
|T1[ □ 



Remark 6.11. One of the important features of the constants C(e, n) (and hence of C^) is 
that they have modulus one. It might be interesting to see that, without any explicit calculation, 
this comes as a consequence of the unitarity of the Fourier transform with respect to the 
norm. As this does not appear - to my knowledge - in the literature, here is the argument : 

Acting on functions of the Schwartz space 5(M^), J- commutes with the operator j-^, and 
hence leaves stable its eigenspaces £n, associated to the eigenvalues 27rn. These are of course 
the subspaces of the Fourier series decomposition in 6 variable : <S„ = {e*"^/(r)}. Here, / is a 
rapidly decreasing function. Therefore, there exists a transformation .F" on such that 

.F(e— V(r))(p,a) = e— ".F"(/)(p). 

J- being an isometry for the norms L'^{M? ,rdrd9) and L'^ (M? , pdpda) , we deduce that J^^ is an 
isometry from ( [0, oo[, rdr) into L^([0, oo[, pdp). 

On the other hand, let Q be the transformation assigning to a function /(r) € L'^{[0, oo[, rdr) 
the function 

roo 

ie—1 . 



G{f){e) = / r''-'f{r)rdr. 







Lemma 6.13 Q is an isometry from ( [0, oo[, rdr) into L^(M, de). 

Proof . The change of variable x = In r allows us to reduce the problem to the unitarity 
of the Fourier transform : let g{e) = G{f){e). We find g{e) = .F(/(e^)e^)(— e), which gives 
hh = Il/(e")e1|2 = C \fir)?rdr = \\fh. □ 
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We know from the definition of C(e, k) that for any function of the form 

<^ = ^(e— V(r)), 
we have : = C{e,k){u(^_^^i^^,if), that is : 

{e''^'r'^-\j^^) = C{e,k){e''^y-^'-\^), 

or : 

Both sides of the equality vanish if n — A; 7^ 0. So let us assume k = n. We obtain : 

gif){e) = i-irC{e,n)gTVi-e), 
which is best depicted by the following commutative diagram : 

L2 ( [0, 00 [,rdr) L^{R,e) 



(-l)"C(£,ri) 



L2([0,oo[,p(ip) 



All transformation involved {Q, T"^ and e — > —e) being unitary ones, one can indeed deduce 
that e I— >■ C(e, n) is a function whose modulus is 1. □ 

A 



Remark 6.12. The above introduced transformations Q and T"^ are related to Mellin's and 
Hankel's transforms M. and Hn respectively (see [14| for their definitions). More precisely, we 
have : 

{Qf){e) = {Mf){s = ie + l) and = r^'Tin. 

(see for instance The previous remark thus yields the following remarkable relation between 

these transformations : 



V/, Mf{i£ + 1) = eC{e, n).{Mnnf){-ie + 1). 



(20) 

A 



Remark 6.13. Formula ( [20| ) allows us to give another proof of the expression of C(e,n). 
First of all, since TL-n = (— l)"Wn, it readily gives C{£,n) = C{£,—n). Now for n ^ choose 

-rl 

fn = r"e 2 . /n is fixed by Hankel's transform of order n. On the other hand, it is easy to see 
that Mfn{s) = 2^-ir(^). This yields : 



Cie,n) 



2^/2r((s + n)/2) 
which is easily identified with formula ([T^). 



, with s = ie + 1, s = —ie + 1, 



□ 
A 
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6.7.3 The outer holonomy 

Let us turn now to the study of the "outer holonomy", that is, the product 

Recah that its quotient by the local holonomy defined in the previous paragraph gives the 
value of the Cech cocycle on . 

Proposition 6.14 • KC^)""'] = 1 + 0{h°°); 

• arg((C'^)°"') € T admits an asymptotic expansion in ^ —1 powers of h, whose terms are 
smooth functions of E; 

• for non zero E, the first two terms in that expansion are the following : 

j ^ «f - ^to + - 4fl argc^ + 

+ ^^K^ + ef;iZn|c| -efiargc + /i(7f)|j , 

where c = efg + iefg . 

Proof. For non zero E, Aq is a regular Lagrangean torus, and we know from proposition 5.6 
that |c^(7f )| ~ 1. This, together with the fact that |C^| ~ 1 too, implies that |(C^)°"*| ~ 1 as 
well. Then of course, since we already knew that (C^)°"* is a smooth function of E, the results 
also holds if ii^ = 0. 

The second point in the proposition is essentially due to our choice of the microlocal solution 
{u^Y on Q£. The distribution 



1 
h 



eh 



{{{x,y),(i,f)))-ef Inp+ef q) I'^C A dr]\ 



2Trih J p 

(by (p, a) we always denote polar coordinates for rj) in the fiber of T*M^) is microlocally on 
equal to a classical oscillatory integral on Aq' H 0^, associated with the canonical Liouville 
half density on Ag', and smoothly depending on E. 



Remark 6.14. This is of course the pendant of remark B.C. Here as well, the phase function 
ip^ = {{x,y), {^,r])) — efolnp + efgO^ is best expressed in terms of complex coordinates zi and 
Z2 '■ we have 

ip^{zi,Z2) = yiiziZ2 - ef lnz2). 

A 

Now, since {u^Y was defined to be the image of that distribution by an classical Fourier 
integral operator, it is also a classical Lagrangean distribution, as are all the other solutions 
{u^y, i 0. We deduce from proposition 5^ and paragraph 5^ that for each i = I, ...,£ — 1, 
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the argument of (c^)*'*"*"^ admits a semi-classical expansion in ^ —1 powers of h whose terms 
smoothly depend on E, as functions with values in T. Hence so does their product (C^)°"*. 

We turn now to the last point of the proposition, which requires a closer look at the principal 
phases of the solutions near the critical point m. Indeed, we know from the non-singular analysis 
of section ^ that the principal phases together with a section of the Keller-Maslov line bundle 
form an invariantly defined object that gives the first two terms of the Cech cocycle X^{h). 

First, let us extend the open sets f^i and to be both equal to \ {m}. This is possible as 
the microlocal solutions {u^Y {u^Y are actually well-defined on all Q. This ensures that 
for every non zero E, IJi=i ^ f^^l cover of the image of 7;p. Still assuming that E is non 
zero, we get from theorem ^]8| that arg(c^(7f')) is, up to 0{h) (and modulo 27r), equal to 



On the other hand, c^iji) is by definition the value of the product 

The problem is thus reduced to the computation of arg(c^)^'^. It is unchanged by conjugation 
by a Fourier integral operator on J7. 

Let us denote by ipf{zi) and ipf{zi, z-i) the phase functions of the oscillatory functions u^e 
and !F^^u_^E ^e . They have the following expressions (modulo 2iThZ) : 

ipf{zi) = <K(ef Inzi) and ipf{zi,Z2) = m{ziZ2 - ef lnz2)- 

As we already saw, for non zero E, they represent - of course - the same Lagrangean manifold 
Aq' = {ziZ2 = Cq'} (the terms fH(eQ'lnii) and 9^(eQ'ln2;2) are actually Legendre transforms of 
each other). Their difference is a function of c = : 

vfizi,Z2)-vfizi)=mic)-^i€^lnc), 

whose real Hessian is the 2x2 matrix whose lines are respectively identified with /c? and 

— E 

— ie^/c^. Its determinant is — If 7^ then eg' 7^ so the signature of that matrix is 

zero. This means that the value of Maslov's cocycle between u,e ,e and J-7^u ,e ,e is zero. 

h ' h h ' h 

The principal phases of those oscillatory functions are expressed by a formula similar to the 
one expressing the phase functions, but with a contribution of the sub-principal terms of 
included. The two principal terms in the argument of (c^)^'^ is their difference : 

arg(c^/'^ ^ I {^(ei) -mie^lne^)) -mieflne^) + 0''ih). 

Here, we have fixed E ^ 0, so the 0^{h) depends on E. 

Finally, together with equation (^l|), we get that the first two terms in the expansion of 
((7^)°"' are, for non zero E : 
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and the proof of proposition |6.14 is completed. □ 



Of course, this also terminates the proof of theorem 3.9. □ 



6.7.4 Recover the regular conditions 

The aim of this paragraph is to show how the quantization conditions of theorem |6.9| , when 
restricted to a non-singular Lagrangean, are equivalent to the regular quantization conditions 
of theorem |5i 



Let us fix 7^ 0. We get a smooth holonomy 1-form [A^] whose action on a basis of the Hi 
of the torus Aq' gives the two quantization conditions. Let 7f' and 7!' be, as previously, loops 
on A^ representing such a basis. 

The loop can be assumed to lie in the neighborhood Q where the normal form applies. 
In the linearized coordinates, one can choose 7^' to be a simple orbit of the X2 vector field. On 
each point of that orbit we saw that the solutions were spanned by 



1 

r 



Therefore, the value of [A^] on 7I' is exactly the monodromy coefficient ^^^^ W ^ quantization 



condition given by theorem 5.5 is thus 



e 27rZ + 0(/i~), 

which is indeed the same as the first singular quantization condition of theorem 

Remark 6.15. This simple observation gives rise, for non zero E, to an intrinsic definition of 
the semi-classical expansion e2{h) introduced in lemma |6l5| , namely : 

4{h) = ^ / A^(M. 
2vr J^E 

Since e|' is smooth ai E = 0, this implies that it is intrinsically defined by the microlocal system 
on a neighborhood of m, that is, independently of the Fourier integral operator U. This will be 
restated in the next paragraph. A 

Let us turn now to the second quantization condition. We know that the action of [A^] on -yf 
is given by the argument of the quotient ((7^)°"*/^'^. Using Stirling's formula for arg(r(z)) = 
3{lnT{z)), with z = -|- ^ (e(/i) = ef{h) + i€2{h); recall that if -E / 0, the principal term cq 
of e(/i) is non-zero), one can perform a very nice cancellation of undesirable terms in the local 
holonomy . 
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Indeed, one finds that 

2 arg(r(z)) = ^ (fH(e In e) - 9^(e) + j(e)^ - ^{e ln(2/i)) ) + 0(/i). 
Here of course, 0{h) depends on E. This gives 

arg(C^) ^ i (9^(elne) - 9^(e)) + 0(/i), 



that is, 



arg(C^) = iin(eolneo -eo) +«n(eilneo) +0(/i). 



We thus find that, if ^ / 0, 

J^E h J^E J^E 2 



which gives the regular quantization condition of theorem 5.8 



7 Structure of the joint spectrum 

This section is devoted to the application of our microlocal study to the precise study of the 
joint spectrum of 2-degree of freedom quantum integrable systems around the critical value of 
a focus-focus singularity. 

More precisely, we shall assume in all this section that we are given two commuting essen- 
tially self-adjoint /i-pseudo-differential operators of order zero Pi{h), P2{h) on a 2-dimensional 
manifold X, with the following properties : 

• the joint principal symbol p = (pi,P2) is proper with connected leaves; 

• p admits a critical point m of focus-focus type, with critical value G M^, such that the 
critical set Aq = p~^{0) has m as its only critical point; 

• at least one of the Pj(/i)'s is a classical elliptic pseudo-differential operator in the sense of 
Hormander (for instance, Pi{h) is a Schrodinger operator : 

2 

for some Riemannian metric g on X). 



Pi{h) = --A,j + Vix), 



The last point implies that the solutions Uh of Pj{h)uh = (j = 1,2) cannot have a wave 
front set at infinity (see [|lO|), which allows us to microlocalize Uh on the Lagrangean level set 
Ao=p-HO). 

Note also that Pi{h) and P2{h) really commute (and not only modulo 0{h°^)). This assump- 
tion is perhaps unnecessary, but we take it in order to use the definition of the joint spectrum 
from ||5j. 

A typical example of this situation would be a Schrodinger operator on a surface of revolution 



with a radial potential, like the quantum spherical pendulum [15| or the quantum Champagne 
bottle 0. As a matter of fact, we will illustrate our results with the latter example, thanks to 
numerics kindly provided by M. Child. 
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7.1 The joint spectrum 



Let be the joint spectrum of the operators Pi{h) and P2{h). We know from Q that the 
assumption of properness for the momentum map p imphes that the intersection of with 
any compact K C M" is a discrete spectrum : each E{h) = {Ei(h), E2{h)) G r\K \s isolated 
and the joint eigenspace 

FE{h) = Mh) G L\X), {Pi{h) - Ei{h)Mh) = {P2{h) - E2{h)Mh) = 0} 

has finite dimension. 



Theorem 7.1 There is a compact neighborhood K of € R" and an Hq > such that all the 
joint eigenvalues E{h) £ n K for h < Hq are simple. Furthermore, the distance between 
two distinct joint eigenvalues is bounded below by some finite power of h. 



Proof. This is a consequence of proposition 6A, apphed to the operators Pf{h) "== Pj{h) — Ej 



Indeed, let X CC M" be a neighborhood of zero such that this proposition holds uniformly for 
E = (Ei,E2) G K. Let ipi{h) and ip2ih) be two normalized joint eigenfunctions for the joint 
eigenvalue E(h), such that (V'l, V'2) = 0. If we can find such eigenfunctions for a set of values of 



h accumulating at zero, then ipi{h) and ip2{h) are both microlocal solutions (see section L5) to 
the equation 

{Pj{h) - Ej{h))i^k{h) ^ i = l,2 A; = 1,2. 

Therefore, their wave front sets are included in a compact W C T*X. Let vr : T*X X he the 
standard projection and let U be an open neighborhood of tt{W) such that U is compact. Then, 

WFh{A)n{T*{X\U)) = (l> (A; = 1,2), 

and because tpk cannot have any wave front set at infinity, -0^ is actually 0{h°°) on X \ U. 
Consequently if we let x be a smooth function with support in U and equal to 1 on tt{W) then 

MkhHx) = m\LHx)+0{h^) = l + 0(/i~). 



and (xV'i,xV'2) = Oih" 
That allows us to replace, in the following, ipk by X'^fc- 



Because of proposition OA, there exists a constant C{h) e Ch such that 

ipi{h) ~ C{h)i^2{h) on T*U. (22) 



Combined as before with the fact that nothing happens at infinity, equation ( |22[ ) implies, by 
taking scalar products, that C{h) simultaneously satisfies |C(/i)| = 1 + 0{h°°) and C{h) = 
+ 0{h°^), which is impossible. This proves that for small h all joint eigenvalues E(h) must be 
simple. 

Now, if E{h) and E'(h) are distinct joint eigenvalues but differing only by 0{h^), the same 
argument applies to the eigenfunctions tpi{h) and ip2{h) associated to E(h) and E'{h), and thus 
proves the last statement of the theorem. □ 
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This theorem extends the results of where it is shown that the eigenvalues are simple 
as long as K does not meet any critical value of the momentum map p. Note that this latter 



statement comes as a consequence of our proposition 5.1 as well. 



Finally, remark that the proof of theorem 7A actually shows that joint eigenfunctions ip{h) 
corresponding to joint eigenvalues E{h) £ K are in one-to-one correspondence with microlocal 
solutions of the system 

{Pj{h) - E,{h))uh ^ Q 

in a neighborhood of UsG-ftT^^' other words, these eigenfunctions are determined by the 



quantization conditions of theorem 6.9. 



7.2 Quantum monodromy 

Away from the critical value of the momentum map p, the joint spectrum locally looks like a 
lattice hi?. However, simply patching these lattices together around is not possible. This 
phenomenon is called quantum monodromy and is described in details in [40|. Let us take 
advantage of section ^ to give a short description of it. 

If Co is regular value of p, then on a neighborhood B of cq there exists a map f{h) : B ^ 
M? that "asymptotically identifies n B with the straight lattice /iZ^", in the following 

sense : f{h) is an elliptic semi-classical symbol of order zero such that for any family E(h) = 
{Ei{h),E2{h)) eB,heH, 

E{h) G f^B + 0{h°°) if and only if f{E{h)- h) G hZ'^ + 0(/i~). (23) 

Indeed, for c G -B let \c{h) be the semi-classical holonomy 1-form on Ac defined in proposition 



5.5; then because of theorem 5.8, one can choose 



/(c;/i)=(:^/ Uh),^ j Ae(/.)), (24) 

J^^(c) Jr2(c) J 

where (71(c), 72(c)) is any smooth basis of Hi{Ac). 

It turns out that up to conjugation of f{h)/h by an element of the affine group with integer 
coefficients GA{2,Z) such a map f{h) is unique modulo 0{h°°). If (5 is a simple loop through 
Co enclosing the critical value and small enough not to enclose any other critical value of the 
momentum map, we can cover its image by a finite number of open ball Bq, . . . , B^. Assuming 
that Bq = B, we can construct a corresponding map fi{h) on Bi that coincide with /(/i) on 
BqCiBi. Continuing this way, we end up with a map fi^h) on B^, and an element /Li((5) G GA{2, Z) 
such that 

^f,(h)=^,{5)o(^f{h)) + 0{h^). 

fi{6) depends only on the homotopy class of 5 and is called the quantum monodromy of the joint 
spectrum. Because the fj{hys can be chosen as in equation (p^, this fj,{6) is actually dual to 
the classical monodromy of the bundle 7ri(Ac) c : 
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Proposition 7.2 ( pO|] ) /f (71(c), 72(c)) and e are chosen as in proposition 6.S then 

where i is the inclusion GL(2,M) ^ GA(2,M) such that for any M E G-L(2,M), t/ie origin 
G is /ixed 6y t(M). 



Tlie article |40] shows how this quantum monodromy can be told from a picture of the 
spectrum by detecting the change in the lattice structure of the joint spectrum around the 
critical value. One of the ways of doing this is to unwind the spectrum onto R^, in the following 
way. Suppose 5 is actually drawn on the universal cover of the punctured plane : TV^!— M^yjO} 
and let Bj be the open cover of its image such that 7r{Bj) = Bj. Then the maps fj{h) on Bj 
define a global map f{h) on i?o U -Bi • • • U B^. By definition we have 

f{s{l)■,h)/h = ^li^)of{s{o)■h)/h. 

Since the restriction of f{h) to the points belonging to the spectrum -7r~^(S(/i)) takes values in 
hZ^ + 0(/i°°), there is for small h a unique map 

f{h) : 7r-\j:{h)) n[jBj ^ hi? 

3 

that coincides with f{h) modulo 0{h°°). This map is of course as well defined on the intersection 
with 7r^^(S(/i)) of any domain T) oiW that projects into a small (but fixed) annulus T) around 
the critical value 0, and is called the unwinding (or developing) map of the spectrum. 

For instance, if one uses this map to unwind a sequence of points Aj{h) G 'k~^ {Ti{h))r\B j with 
TT{Ai{h)) = 7r(Ao(/i)) (see figure P), then one obtains an integral open polygonal line (Aq, . . . , Ai) 
in hi?, whose extremities satisfy 

A,/h = ii{5){AQ/h). 

For short, we shall allow us to refer to such a sequence of points Aj as a closed polygonal line in 
S(/i) whenever these points are close enough to each other so that the line segments [^j,74j-|_i] 
lie in /(P; h). The "segment" between two consecutive points Ajiji) and Ajj^iiji) is by definition 
the preimage in V of [74^, 



7.3 The exact counting function 

The purpose of this paragraph is to count the number of joint eigenvalues inside a simple 
closed polygonal line around the critical value 0. Let Lo(^) be the subset of S(/i) n I? of joint 
eigenvalues E{h) that are unwound on the horizontal line {y = 0} C (in other words, 
f2{'^~^{E{h));h) = 0). Because of proposition |7.2| , LQ{h) is pointwise fixed by the quantum 
monodromy fj,{6). That means that any closed polygonal line A{h) = {AQ(h), . . . , Ai{h)) on 
S(/i) starting on -Lo(^) is unwound as a closed polygonal line A{h) = {AQ{h), . . . ,A£{h)) in 
(see figure ^) . 

It is easy to see that if we have another such polygonal line on S(/i) n D that lies inside the 
previous one, then the number of joint eigenvalue between these two is exactly the number of 
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^3 



Figure 5: unwinding of a sequence of points in H V, in the case of the spectrum of the 

Champagne bottle (numerics by M. Child). In the case h), the starting point belongs to LQ{h), which 
implies that the unwound line is closed. 
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^-integral lattice points between the two unwound polygons (see lemma . As a consequence, 
there should be a universal index v ^ TL such that the number A^/j(^(/i)) of joint eigenvalues 
inside A{\i) be 

where Nh.{A.{\i)) denotes the number of /i-integral points inside the polygon AiK). Note that 



the latter is given by Pick's formula |36] : 



N.^Am = A.e.iAm + ^^--^^^^(^(h)) + 1, (25) 

where ZLength(^(/i)) is the number of /i- integral points on the boundary. 
The result is actually that = : 

Theorem 7.3 Let A{h) = (AQ{h), . . . , be a simple closed polygonal line on starting 

on Lo(/i) and enclosing one the origin, and let A{h) = (AQ^h), . . . , Ag{h) = Ao(/i)) be the 
unwound polygonal line in E?. Denote by Nh{A{h)) the number of joint eigenvalues inside A{h) 
(counting the points in the boundary), and similarly denote by Nh(A{h)) the number of h-integral 
points inside the polygon A{h), counting the boundary. 

Then 

Nh{A{h)) = Nh{A{h)). 

In other words, the recipe to calculate Nh{A{h)) is simply to unwind A{h) and then apply Pick's 



formula (25). 



Proof. The idea is to cut the "polytope" delimited by A{h) into two halves separated by Lq; 



then on each part the control over the joint eigenvalues given by theorem 6.9 is used to shift the 



problem into a region where the result will be obvious, due to the following fact : 

Lemma 7.4 if B is a simply connected open subset of T>, then for any closed polygonal line 
C{h) in Ti{h) H B, the unwound line C{h) is also closed and we have 

Nh{C{h)) = Nh{C{h)). 

This lemma also proves the existence of the index i' claimed before the statement of the theorem, 
and this implies that we can choose any particular A{h) to prove the theorem. So let us assume 
that the intersection of the polytope delimited by A{h) with the horizontal lines €2 = hn are 
segments whose extremities are vertices of A. 



Now, let us apply theorem 6^ to the operators Pf{h) =^ Pj{h) — Ej, and keep the notations 



from the statement of this theorem. We know from lemma |6.5| that the map 

E = {E^,E2)^e^ih) = iefih),eiih)) 

is an elliptic semi-classical symbol of order zero. It is therefore for small h a diffeomorphism 
of a compact neighborhood K of zero into its image; this implies that the determination of the 
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joint spectrum in K amounts to that of the corresponding values of e^{h), which we shall hence 
simply denote by e = (ei, £2)- From now on, we shall assume that the annulus P defined in the 
previous paragraph is a subset of K. 

Let us look for solutions satisfying €2 ^ 0. Let gi{h) and g2{h) be the following functions : 

— 51 e;/i = A^r' + e2--i^ln 2/1 -2argr' ^' 
n In 



and fif2(e; h) = €2. 

Because of the analyticity properties of L, these functions are holomorphic in €2^0. The map 
f{E] h) '= g{e^{h); h) satisfies equation (p3|) thanks to theorem |6.9|. Because of the In/i term in 



it, it is not a standard semi-classical symbol globally; but we know from paragraph |6.7.4| that it 



is so for S in a small compact at finite distance from the origin. Moreover, because of remark 



6.15 , it must satisfy equation ( P4|) as well, where (71(c), 72(c)) is the basis given in proposition 



6.3. 



Therefore, the joint spectrum is spread on the horizontal lines €2 = hn, n £ 1^, and the set 
LQ{h) lies in the line 62 = 0. Our problem is now to count the number of joint eigenvalues inside 
the intersection A'^{h) of the half plane €2^0 with the polytope delimited by A{h) (figure 
Let us denote by In{h) the set of joint eigenvalues contained in A'^{h) and satisfying €2 = hn. 




Figure 6: the upper half of the polytope. 



From now on, we forget the negative half of our polytope, and label 

io,...,i£ 

the consecutive vertices of ^+(/i), starting and ending on Lo(/i). We denote by Nh{A{h)) 
the number of joint eigenvalues inside A'^{h), and by Nh{A{h)) the number of /i-integral points 
inside the polytope delimited by the unwound points Aq, . . . , Ag, Aq. The theorem will be proved 
provided we show that Nh{A{h)) = Nh{A{h)). 

The next step of the proof is now to translate A~^{h) far from the origin, as follows (figure 0). 
The domain T> where the unwinding map is defined can always be chosen such that there exists 



60 



an integer k{h) £ TL such that the origin is outside of the polytope C(K) = (Co(/i), . . . ,Ci{h)) 
defined by 

fiCj {h);h) = f{A,{hy,h) + hik{h),o). 




Figure 7: translation of the polytope 



Of course, the set of joint eigenvalues inside C{h) satisfying €2 = hn is in bijection with /n(/i), 
which imphes that 

Nh{C{h)) = NhiMh))- 

On the other hand, because f{-',h) is a vahd unwinding function in a region containing 
both A{h) and C{h), the unwinding C{h) of C{h) is by definition the translation of A{h) by the 
horizontal vector h{k{h),0). It has therefore the same number of interior points : 

NhiC{h)) = NhiAih)). 



But lemma 7.4 yields 



Nh{C{h)) = Nh{C{h)), 



which finally proves the theorem. 



n 



7.4 Shape of the spectrum near the critical value 



We investigate here the shape of the joint spectrum in a region of size 0{h) near the origin, in 
terms of the (ei,e2) variables. We already know that the joint eigenvalues are distributed on 
the horizontal lines e2 = hn, n G Z. It remains to study their repartition on each of these lines. 
For this purpose, we fix n E Z and denote by x = €i/h the generic variable on {e2 = hn}. Of 
course, the results will be based upon corollary 6.10| , whose notations we use here, and which 
we express in the following way : 
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Proposition 7.5 (corollary |6.10 ) Let K a compact neighborhood of the origin in'M? . For any 
family E{h) = {Ei{h), E2{h)) G (HK) such that e^^^^h) = hn + 0{h°°) (n e Z), the following 
characterization holds : 

E{h) G + 0{h°°) if and only if 



|n|^ -x(/i)ln(2/i) -2argr J^!^W^i±H^ + /i) ^ 27rZ + 0(/i= 

with x{h) = ef^'^\h)/h. Here fn{h) is a symbol of the form : 

fn{x; h) = ^ + xB + nC + D + fi^ + /i/n,i(x) + h^fn,2{x) + ■■■ , 
for some real constants A,B,C,D. 



Proof. Of course, in view of corollary 6.10| - applied to the system Pj = Pj — htj, j = 1,2 -, A 
is just the action integral on the critical Lagrangean : A = J^o ao, and it only remains to prove 



that there exists constants B, C, D such that 



1^0 (k*) =x{h)B + nC + D + 0{h), 



where t and x{h) are linked by 



x{h) = e\\h) = ]-ef{h). 



Because the sub-principal symbols of are affine functions of t : f* = r — t, it is easy to deduce 
from formula (|l7|) of proposition that /^o(k*) is affine in t. The result then follows from the 
fact that the principal part of e*(/i) = {x{h), n + 0(/i°°)) is also an - invertible - affine function 
oft. □ 

We study the distribution of joint eigenvalues on the line 62 = hn by investigating the gaps 
between two consecutive eigenvalues. 

Theorem 7.6 For any n G let Ek{h) be the sequence of joint eigenvalues in {e2 = hn}(~\{hK) , 
such that the corresponding variables Xk{h) form (for each h) a strictly increasing finite sequence. 

The gap between Ek{h) and Ek+i{h) is given by : 
Ek+iih) - Ek{h) 



h 



lira 

\lnh\ + B-ln2-^'^x; 



il + Oih)) 



(26) 



27r 



lnh\ ' (ln/i)2 ^^"^-^^ " ' '"-^ J ■ " \{lnhf 

for an x £]xk, Xk+i[, and a = \\M ■ (1, 0)||; we have denoted 

T / \ r, 1'^^ + 1 + |n 
^n{x) = 2argr 
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Note that is related to the logarithmic derivative %[) of the Gamma-function : 

Ki.x) = ^[i^ ^ J] . 

A consequence of this theorem is that the number of Ek{hys inside hK increases as h tends 
to zero at a rate equivalent to const. \ ln/i|, which easily yields a Weyl type asymptotic formula 
for the counting function (see theorem [7.7| in next paragraph). A more precise result is that 
the function controls the shape of the graph "gaps versus x" . It turns out that has a 
non-degenerate minimum at x = 0, which is particularly sharp for n = (figure ^). This might 



*n(^) 




Figure 8: the function ^o(x) (here 7 denotes Euler's constant 7 ~ 0,5772156649 . . . ). 



have some interesting applications in numerical analysis. First, it should refrain people from 
thinking that the gap curve has a cusp; secondly, however, this "near-cusp" might be useful to 
precisely locate the critical value a; = 0. Numerics for the example of the Champagne bottle 
suggest that both remarks apply for reasonable value of the small constant h (figures |9| and 10). 
Finally, of course, this theorem can be used to find a rigorous approximation for the spectral 




-0.2 



Eh 



Figure 9: spectral gaps for the quantum Champagne bottle, n = 0, /i = 10 ^ (numerics). 



gaps. In the example of the Champagne bottle, formula ( p^ ) can be made completely explicit. 
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The Hamiltonian H is the following : 

H{x, y, r/) = ^(e^ + r/2) - + r^ 

with = + y^. It commutes with the angular momentum / = xrj — y^, hence we can let 
p = {pi,P2) = {H,I). The corresponding /i-operators are 

A 2 4 t h d 
2 z 



One can compute the quantities involved by formula (26), and one finds 

y/2 



" ■ 



J ^ SO o = \/2, and 5 = (5/2) In 2. 



Thus, using formula (|2q ) to approximate the smallest gap (obtained at x = for n = 0) we get : 

. ( Ek+i-Ek \ 27r^/2 , ^^^^ 

k \ h J |ln/i| + (9/2)ln2+7 ^ ^ 

The comparison between the principal part of this formula as a function of In h and the "exper- 
imental" result from the numerical computation of the spectrum is depicted in figure 11, and 
tends to show that the approximation is quite accurate. 



Proof of theorem ^ Because (Ek+i - Ek)/h = {M + 0(h)) ■ {xk+i - Xfc, 0) + 0(/i~), it 
suffices to estimate \xk+i — Xk\- 

For n G Z, let 

gn{x;h) =^ ^ ^|n|^ - xln{2h) -2argr ^ ^^ + ^+ j _^ ^^(^../j) 

The function x gnix; h) is smooth and we have 

g'^ix; h)=B- \n{2h) - <(x) + 0{h). 
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For small h, gn{h) is therefore strictly increasing, with slope of order | ln/i|. The solutions Xk{h) 
to g{xk{h);h) = k £ + 0{h°°) are therefore, for small h, in one-to-one correspondence with 
the exact solutions to g{xk{h);h) = k (except maybe for the first and the last ones, sitting 
near the boundary of the compact interval in which we are looking for solutions). Moreover 
the differences between the exact and the 0{h°°) solutions are of order which is a 

fortiori 0{h^). In the following denotes the exact solution. 

By RoUe's theorem, there is an x £\xk-,Xk+i[ such that 

i =g'^{x;h). 

\Xk+l - Xk\ 

This means 



\lnh\ +B-ln2- ^'^(x) + 0{h) ' 
which proves the theorem. □ 



7.5 Weyl's formula 

Let K he a. compact convex subset of R^, and let K{h) be its image by the change of variables: 

Recall that the principal part of this map is 

el = ■ {t - r(m)), 

where r(m) = (ri(m), r2(m)) is the value at m of the joint sub-principal symbol of the system 
Pi{h), P2{h), and M = dF{0), where {pi,P2) = FiQi,Q2) for some symplectic coordinates near 
m (see section 6.5). 
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Denote by Kq the principal part of K{h); it is a convex compact of M^. Theorem 7.6 easily 
yields the following estimate for the counting function : 

Theorem 7.7 Denote by Nf^(K) the number of eigenvalues belonging to hK : 

NhiK)''^' #{{E,,E2) e hK n^{h)}. 

The following estimate holds, as h tends to zero : 

\nh\ 



Nh{K) 



2tt 



\dx ^dn\ +0(1). 



The measure \dn\ is the counting measure on Z, and \dx\ is the standard one- dimensional 
Lebesgue measure. 

We present in figure ll^ a numerical illustration of this result. 




In h 



Figure 12: The counting function Nh{K) for the Champagne bottle, where is a compact selecting 
only joint eigenvalues on the line 62 = (numerics). 



The question that this theorem arouses is whether the formula for A'^^ can be related to a 



Weyl-type heuristic involving volume in phase space. The answer is positive (proposition 7 



up to the difference between the measures of Kq with respect to \dx f\ dn\ and with respect to 
the standard Lebesgue measure on R^, which of course strongly depends on the geometry of Kq. 

Let IX be the push-forward on of the symplectic (or Liouville) measure by the momentum 
map p = {pi,P2) '■ for any Borel B of M^, 



KB) 



Iw Aw|/2. 



Let \dc\ = \dci A dc2\ be the standard Lebesgue measure on M?. The Lebesgue volume of B will 
be denoted by \B\. 
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Proposition 7.8 As h tends to zero, 

Proof. Let VLhe a, neighborhood of the critical point m in which the system (^1,^2) can be 
brought to the normal form of lemma 6^, so that there is a local diffeomorphism F defined on 
U = piyt.) such that, with the same notations as in section 6.2, (^1,^2) = F{qi,q2). 

Replacing K by hoK with ho small enough, if necessary, we can always assume that K C U. 
Let e > be such that the ball 

Be'^= {{r,9,p,a) £ M, r<e,p<£} 

is included in Q. Then iJi{hK) can be split into two parts : 

fi{hK)= I \dxdyd^dr]\+ [ |a;Aa;|/2. (27) 

JB,n{iqi,q2)<^F~HhK)} J B^^n{ipi ,P2 )ehK} 

If h is small enough, each fiber p~^{c), c G hK intersects B^, and we know from section 
|6^ that the fibration p^B^ over hK is trivial with fibers diffeomorphic to compact cylinders. 
Moreover, p naturally defines on each cylinder a locally invariant measure p such that the 
symplectic measure splits according to : 

\ll> A oj\ / 2 = p (g) \dpi A dp2\. 

Therefore, the second term of (^) can be written as : 

/ |a;Acj|/2= / p{B^ np'^{c))\dci Adc2\. 

J Bin{{pi,p2)&hK} J{ci,C2)ehK 

Because p{B^ np^^(c)) is bounded for bounded (ci,C2), the following estimate holds, for some 
C > : 

[ \ujAuj\/2<C. [ \dciAdc2\=C.h^\K\. 

JBin{(pi,p2)<^hK} J{ci,C2)<^hK 

Let us turn now to the first term in the sum (UTi) : 



Ie{h) '= / \dxdyd^dri\. 

JB^nUcn .a2)eF-l(hK)\ 



I B,n{{guq2)eF-HhK)} 

Using the new coordinates {r,6,i,X) with i = rp and X = a — 9, the symplectic measure is 
transformed into 

i 1 

-\dX A di A d9 A dr\ = -\dqi Adq2\ \d9 Adr\. 

This gives : 

Ie{h) = 27r / l|H<e2 I [ , -\dr\ ] \dci A dc2\ 

J{cuC2)&F-HhK) V-^^ J 
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f £2 

2vr / l|c|<e2 ln(— )|cici A (ic2|. 



-2 

Note that for h small enough, the condition |c| ^ is always satisfied for c G F~^{hK). 



Now, the change of variables c = F{c)/h yields 

lAh) = 2^h^ [ H f )J{hc)\dc\, 
JceK 1^ {hc)\ 



where J is the Jacobian of F ^. Taylor's formula for J and F ^ easily shows that the integrand 
can be reduced to its principal part : 

I,{h) = 27rh^ [ H / )jmdc\ + 0{h^). 
JceK \M ^.hc\ 

Since J(0) = detM~^, the right hand-side is equal to : 

27r/i2 / in(^)|dc| +0(/i3), 

with c = M~^.c. Because the measure ln|c||dc| is locally integrable at the origin, this finally 
gives : 

= -27Th'^lnh\M~\K)\ + 0{h^), 
thus proving proposition |7.8| . □ 
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